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IX 


Incroduct ion 


I  . 

Refraction  data  collected  at  the  Lajitas  site  in  1983 
yielded  in  situ  P-wave  velocities  of  3.5  km/s  for  the  Santa 
Elena  limestone.  A  regional  refraction  survey  conducted 
with  a  Vibroseis  and  a  receiver  line  extending  north  from  the 
*  seismic  station  to  just  a  few  miles  south  of  Alpine  alo.ng 

Highway  118  yielded  P-wave  velocities  of  3.5  km/s  for  the 
first  160  meters  below  the  earth's  surface  and  4.7  km/s  for 
the  next  1350  meters  (Figure  1.1)  (Golden  et  al.,  1  985)  . 

Three  local  refraction  surveys:  1)  an  1100  .meter  long 
reversed  refraction  spread  with  a  dynamite  source  (line  ori¬ 
ented  NNW/SSE) ,  2)  a  38  meter  long  reversed  refraction  spread 
with  a  sledge  hammer  and  aluminum  plate  source  (line  oriented 
E/W)  ,  and  3)  a  down-hole  survey  at  a  IGO  m.eter  depth  with 
sledge  ham.mer  and  plate  source  (shots  at  three  meter  inter¬ 
vals  extending  out  to  approi-'.imate ly  80  m>'‘,:ers  south  of  hcre- 
hole)  yielded  average  P-wave  velocitie.s  of  3.4i’7  ±  C.G5  k:"  s , 
^  3.2  +  0.05  km/s,  and  a  mean  apparent  P-wave  velocity  of  3.4"^" 

km/s  (deviation  =  ±0.16  km/s)  respect  iv-v  ly  (Figur':'  '..2) 
(Reinke  and  Logan,  1983) . 

Based  on  the  data  accumulate^!  hy  Ferre)  )s  (1 
den^r  (’91’4),  an-i  Clark  (1  966)  cer.c  ' r.  i  re;  P-ven'”  r'- hr  :  :  ‘  e  r 
in  hom-rafr  r:eoi_e;  1  i  mrcrt ,  the  in  si'n  F-wavr  ■'Ur- hre  : r 


2 


3 


re  1.2  Location  of  local  refraction  survey  traverses.  North  points  to 
top  of  the  page  and  1  cm  equals  approximately  0.23  km.  {From  Reinke  an^ 


the  Santa  Elena  limestone  were  expected  to  average  approxi¬ 
mately  5.5  km/s  (Figure  1.3). 

Sorrells  measured  P-wave  velocities  through  three  mutu¬ 
ally  perpendicular  cylindrical  plugs  cut  from  each  of  four¬ 
teen  carbonate  rock  samples  from  different  localities.  Each 
set  of  plugs  measured  2.54  -  7.62  cm  in  length  and  2.54  cm  in 
diameter.  For  wavelengths  between  10  and  20  millimeters,  the 
measured  longitudinal  wave  velocities  through  each  plug 
varied  with  changing  hydrostatic  pressures  ranging  from  1  - 
2000  bars.  At  one  atmosphere,  velocities  varied  from  4.65  to 
6.36  km/s  in  the  samples  and  as  much  as  thirty  percent 
azimuthally  in  individual  samples.  Sorrells  also  examined 
the  variations  in  the  petrographic,  physical,  and  min- 
eralogical  properties  of  the  fourteen  rock  samples.  Based  on 
correlations  between  the  velocity  data  and  the  petrographic, 
physical,  and  mineralogical  variations  Sorrells  (1961) 
concluded  that:  1)  relatively  large  decreases  in  velocity 
correlate  with  relatively  minor  increases  in  porosity  ranging 
from  0-3  percent;  2)  increased  pressure  reduces  the  poros¬ 
ity  thus  causing  a  velocity  change;  3)  a  range  of  0  -  30  wt ' % 
in  clay  content  decreases  P-wave  velocities;  4)  an  increase 
in  dolomite  content  produces  a  slight  increase  in  longitudi¬ 
nal  wave  velocity;  5)  linear  rock  fabrics  visible  in  thin 
section  possess  anisotropic  properties;  and  6)  bulk  density 
increases  as  P-wave  velocity  increases . 
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Excluding  samples  with  high  clay  content  and  vuggy 
porosity  yields  a  bulk  density  range  of  2.60  -  2.68  gm/cc  and 
corresponding  P-wave  velocities  ranging  from  5.0  -  6.5  km/s 
for  the  Sorrells  data.  Gardener  (1974)  documents  that  in 
situ  velocities  increase  with  increasing  density  and  in  the 
2.60  -  2.68  gm/cc  bulk  density  range  in  situ  P-wave  velocity 
ranges  from  approximately  18000  -  21000  ft/s  (or  5.49  -  6.40 
km/s) .  Clark  (1966)  also  shows  that  in  situ  P-wave  veloci¬ 
ties  range  from  5.5  -  6.5  km/s  for  homogeneous  limestone. 
Therefore,  in  situ  velocity  measurements  correspond  to  labo¬ 
ratory  velocities  for  uniform  homogeneous  limestones. 

Laboratory  velocity  measurements  vary  significantly  with 
the  nature  of  saturating  fluids  (Gardener  et  al.,  1974;  King, 
1966;  Wyllie  et  al.,  1956,  1958;  Elliot  and  Wiley,  1975);  and 
because  laboratory  observations  reveal  that  P-  and  S-wave 
velocities  depend  strongly  on  porosity  and  saturation  condi¬ 
tions,  Toksoz  et  al .  (1976)  raise  several  questions; 

1)  What  effects  do  pore  geometry  and  saturating 
fluids  (such  as  gas,  water,  or  oil)  have  on 
seismic  velocities? 

2)  Can  the  pore  geometry  and  saturation  state  be 
determined  given  the  seismic  velocities? 

3)  Given  these  rock  models,  for  seismic  velocities 
under  reservoir  conditions,  how  might  one 
determine  from  seismic  data  the  nature  of  the 
saturating  fluids  (e.g.,  gas,  oil,  brine)? 

Several  theoretical  models  have  been  formulated  to 
answer  these  questions.  Biot  (1941)  introduced  a  semi- 
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intuitive  formulation  of  the  equations  of  elasticity  for  a 
porous  aggregate.  Most  importantly,  he  recognized  that  in  a 
solid,  the  pore  fluid  pre.S3ure  and  volume  increment  are  state 
variables,  in  addition  to  stress  and  strain.  Gassmann  (1951) 
expressed  these  variables  in  terms  of  the  separate  properties 
of  the  pore  fluid  and  solid  material.  In  1981,  Burridge  and 

•  Keller  carefully  derived  the  basic  equations  first  assumed  by 

Biot  and  Gassmann  and  clarified  the  corresponding  valid 
frequency  regimes.  A  number  of  models  (e.g.,  Eshelby,  1957; 
Bristow,  1960;  Kuster  and  Toksoz,  1974;  O'Connell  and 
Budiansky,  1974;  Budiansky  and  O'Connell,  1976)  sharing  the 
minimal  assumptions  of  Biot-Gassmann  theory  concerning  pore 
geometry  emerged,  attempting  to  provide  stronger  predictions 
on  the  elasticity  of  aggregates  through  stronger  assumptions 
about  the  microscopic  geometry  of  the  constituents.  Thomsen 
(1985)  extended  the  standard  .Bu  ;  ansky  and  O'Connell  (1980) 
model  theory  for  dii.ute  concentrations  of  fluid 
heterogeneities,  to  high  concentrations  of  pores  and 

^  fractures,  and  simultaneously  preserved  explicit  consistency 

with  the  predictions  of  B lot -Gassma.nn  theory  for  low- 

% 

frequency  elastic  moduli  of  porous  rocks.  The  predictions  of 
Biot-Gassmann  theory  are  as  follows;  1)  the  shear  modulus  of 
an  unsaturated  rock  equal.s  that  of  the  same  rock  saturated 
with  liquid,  and  2)  the  difference  between  the  unsaturated 
and  saturated  bulk  modulus  is  a  defined  amount.  Laboratory 
observations  of  the  "Biot  wave"  provided  support  for  the 
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validity  of  other  predictions  made  by  Biot-Gassmann  theory 
(Fiona,  1980) .  These  other  predictions  relate  unsaturated 
moduli  to  saturated.  Formal  consideration  of  these 
theoretical  models  will  be  delayed  until  chapter  V. 

Fractures  in  the  Santa  Elena  limestone  lower  the 
expected  wave  velocities  in  the  first  160  meters.  The  objec¬ 
tive  of  this  research  is  to  investigate  the  effects  of  frac¬ 
tures  on  seismic  wave  velocities  at  the  Lajitas  site  with 
particular  emphasis  on  macrofracture  density  and  orientation 
and  their  effects  on  in  situ  P-  and  S-wave  velocities.  In 
determining  these  effects,  we  examine  the  relationship 
between  several  measurements  of  fracture  orientation, 
porosity,  and  density  and  in  situ  P-,  Sv-,  and  Sn-wave 
velocities  at  varying  azimuths.  Matrix  and  whole  rock 
parameters,  estimated  from  laboratory  measurements  of  P-  and 
S-  wave  velocities,  pore  porosity,  dry  bulk  and  saturated 
rock  densities,  and  grain  density  for  the  Santa  Elena 
limestone,  aid  in  examining  this  relationship.  The 
laboratory  and  in  situ  velocity  measurements  provide  the 
parameters  necessary  for  the  application  of  the  theoretical 
models  relating  pore  porosity,  fracture  porosity  and  fracture 
density  to  seismic  wave  velocity. 

First,  a  brief  description  of  the  geologic  setting  is 
given.  Then,  a  discussion  concerning  the  spatial  distribu¬ 
tion  of  fractures  at  the  macroscopic  scale,  with  emphasis  on 
macrofracture  density  and  orientation  follows  in  section 
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three.  Next,  attention  is  focused  on  microscopic  inhomo¬ 
geneities  (such  as  vuggy  porosity,  preferentially  oriented 
microfractures,  compositional  variations,  etc.)  possibly 
affecting  in  situ  seismic  velocities.  The  term  macroscopic 
refers  no  structure  and  inhomogeneities  easily  observed  in 
outcrop  with  the  unaided  eye,  to  regional  structure  and 
lithologic  changes  spanning  terrain  distinguishable  from 
aerial  and  satellite  photos  with  scales  averaging 
approximately  1:22,000.  Structural  and  mineralogical  hetero¬ 
geneities,  seen  only  through  a  microscope,  to  hairline 
fractures  and  grain-sized  mineralogical  components  barely 
visible  to  the  unaided  eye,  define  the  term  microscopic . 
Section  five  investigates  the  relationship  between  in  situ  P- 
and  S-wave  velocity  and  fracture  orientation  and  utilizes  the 
Biot-Consistent  model  to  estimate  the  average  fracture 
porosity  and  density  along  each  seismic  refraction  spread  in 
the  survey  area.  Finally,  section  six  contains  a  brief 
summary  and  conclusions.  Appendix  A  shows  the  association 
between  all  figures  containing  maps. 
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II.  Geologic  Setting 


The  region  of  concern  encompasses  the  Lajitas  seismic 
station  located  approximately  32  km  northeast  of  Lajitas,  a 
small  village  near  the  Rio  Grande  river  in  southwest  Texas 
(Figure  2.1) .  Situated  in  the  northeast  quarter  of  the  7  1/2 
minute  Amarilla  Mountain  Quadrangle,  the  study  area  covers 
approximately  5  km^ .  The  region's  southern  border  begins 
about  1  km  north  from  FM  170.  The  edge  of  the  Long  Draw's 
western  floodplain  marks  the  region's  eastern  border. 

The  regional  structure  can  best  be  described  as  a  frac¬ 
tured  antiform  caused  by  a  lacolithic  intrusion  of  unknown 
depth  during  the  Tertiary  (Herrin,  personal  communication) . 

Cn  a  macroscopic  scale  the  predominant  exposed  rock  in  the 
region  is  lithologically  homogeneous  Santa  Elena  limestone 
(Figure  2.2) .  Examination  of  infrared  satellite  and  air 
photos  reveals  a  few  grabens  and  a  synclinal  structure  with 
most  recent  units  composed  of  macroscopically  unfractured 
shale.  Along  the  perimeter  of  the  study  area,  shale  c^estas 
that  tilt  upward  toward  the  central  part  of  the  region 
encircle  the  exposed  limestone  terrain.  The  Amarilla 
Mountain  7  1/2  minute  topographic  quadrangle  shows  a  gradual 
rise  in  elevation  where  the  shale  beds  flank  the  Santa  Elena. 
The  change  in  elevation  increases  where  Santa  Elena  limestone 


Figure  2.1  Location  of  study  area 


is  exposed.  Just  outside  the  immediate  vicinity,  a  few 
Tertiary  intrusions  puncture  the  horizontal  limestone  beds  at 
the  surface.  Lacolith  intrusions  upwarped  limestone  at  Black 
Mesa  (Figure  2.2)  and  probably  caused  the  Soltario  uplift 
located  about  14  miles  to  the  northwest  (Maxwell,  1971) . 

Fractures  cut  the  Santa  Elena  Formation  and  possibly  the 
Sue  Peaks  and  Del  Carmen  Formations  located  strat igraphically 
below  the  Santa  Elena  Formation.  The  air  and  infrared  satel¬ 
lite  photos  re/eal  d.stinct  fracture  systems  of  undetermined 
depth  riddling  the  exposed  Santa  Elena  limestone  (Figure 

2.3)  .  Th^  fractures  could  be  open  to  the  water  table  about 
1500  feet  below  the  surface.  Because  of  the  region's  hot  and 
and  climate,  it  is  likely  that  these  m.acrofractures  are  air- 
filled. 

Shale,  marl,  and  thin  marly  limestone  ledges  make  up  the 
Sue  Peaks  Formation  a.na  massive  L.hick  bedded  limestones 
.mainly  compose  the  Santa  Elena  and  Del  Carmien  Formiations 
(.Maxwell,  1971)  .  Local  refraction  surveys  conducted  in  the 
Lantas  area  show  that  the  Santa  Ele.na/Sue  Peaks  a.nd  Sue 
Peaks  Pel  Carm.en  interfaces  reach  depths  of  approximately  245 
and  320  -  33  5  m.eters  below  the  surface  respectively  (Figure 

2.4)  (Alden,  1983;  Reinke  and  Logan,  1983)  .  All  three  'orma- 
ticnc  are  located  within  the  Com.anche  Series  (Maxwell,  1971)  . 
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P-  &  S-WAVE  VELOCITY  STRUCTURE 
LAJITAS,  TEXAS 
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Figure  2.4.  The  approximate  stratigraphic  and  seismic 
profile  at  the  Lajitas  site  in  southwest  Texas.  (Drawn 
from,  Alden  (1983)  .) 
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Ill . 


Preliminary  Map  and  Field  Work 


The  main  preliminary  map  and  field  work  objectives  are 
to  1)  determine  fracture  distribution  and  orientation  both  at 
a  macroscopic  and  microscopic  scale;  and  2)  collect  oriented 
Santa  Elena  hand  specimens.  Knowledge  of  fracture  orienta¬ 
tion  and  distribution  helps  in  defining  possible  anisotropic 
symmetries  and  consequently  aids  in  planning  refraction  sur¬ 
veys  such  that  proposed  symmetries  can  be  verified  or  dis¬ 
carded.  Comparison  of  laboratory  measurements  of  physical 
properties  in  oriented  hand  specimens  aids  in  defining  the 
homogeneity  and  elasticity  of  the  Santa  Elena  rock  matrix. 

Measurements  of  fracture  orientation  at  both  a  micro¬ 
scopic  and  macroscopic  scale  aid  in  estimating  possible 
anisotropic  symmetries  defined  by  open  fracture  sets.  An 
isotropic  and  homogeneous  rock  matrix  possessing  open  and 
aligned  vertical  fractures  restrains  effective  seismic 
anisotropy  to  a  simple  hexagonal  symmetry  (Crampin  et  al., 
1984;  Crampin,  1984) .  With  the  axis  of  symmetry  oriented  at 
right  angles  to  the  vertical,  this  model  can  be  thought  of  as 
a  transversely  isotropic  medium  rotated  90  degrees  from  the 
horizontal.  Note  also  that  subvertical  fractures  cause  neg¬ 
ligible  change  in  the  symmetry  just  mentioned  (Booth  et  al., 
1986)  .  If  anisotropic  symmetries  defined  by  microscopic  a.no 
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macroscopic  fractures  mirror  one  another  as  Booth  et  al . 

(1986)  claim,  then  anisotropic  symmetries  determined  in  the 
lab  accurately  predict  corresponding  anisotropic  symmetries 
in  situ;  provided  the  rock  matrix  is  isotropic  and  homoge¬ 
neous  . 

Testing  the  homogeneity  of  the  Santa  Elena  rock  matrix 
requires  the  collection  of  several  oriented  rock  specimens 
from  random  locations  within  the  vicinity  of  interest.  A 
comparison  of  physical  rock  properties  for  each  collected 
specimen  aids  in  determining  the  'degree*  of  homogeneity. 
Physical  properties,  such  as  porosity,  density,  mineralogical 
composition,  and  P-  and  S-wave  velocity  can  be  measured  in 
the  laboratory.  Measuring  P-  and  S-wave  velocity  for  each 
sample  along  several  propagation  paths  differing  in  orienta¬ 
tion  aids  in  determining  how  the  physical  properties  affect 
the  composite  elastic  constants  for  the  medium. 

Planning  the  geometrical  arrangement  of  seismic  refrac¬ 
tion  lines  requires  the  knowledge  of  tra^'^ure  distribution 
and  orientation.  To  observe  slight  anisotropic  effects  in 
isotropic  and  homogeneous  material  (i.e.,  a  2  percent  differ¬ 
ence  in  wave  velocity)  ,  a  source  must  generate  wavele.ngt.hs 
much  greater  than  the  fracture  size  and  the  separation  dis¬ 
tance  between  fractures  (Backus,  1962)  .  Because  longitudinal 
wave  velocities  decrease  to  a  minimum  as  the  angle  of  inci¬ 
dence  to  a  fracture  plane  approaches  90  degrees,  one  must 
consider  fracture  orientation  when  planning  the  geometrical 
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arrangement  of  seismic  refraction  lines  (Lynn  and  Thomsen, 
1986) .  Also,  because  tensional  and  compressional  stresses 
open  and  close  fractures,  aligning  refraction  lines  parallel 
and  perpendicular  to  traces  of  vertical  fracture  sets  pro¬ 
vides  a  means  to  determine  the  current  stress  regime  (Rai  and 
Hanson,  1986) .  Knowledge  of  accessible  areas  that  minimize 
terrain  corrections  and  geophysical  field  worlc  is  also 
valuable  . 


Methodology 

Determination  of  Macrofracture  Distribution  and  Orienta¬ 
tion.  Compilation  of  a  detailed  lineation  map  representing 
fracture  lineations  at  the  surface  by  inspecting  stereo  pairs 
and  a  high  resolution  infrared  photo  enables  the  measurement 
of  macrofracture  distribution  and  orientation.  Observations 
in  the  field  both  above  and  below  the  surface  also  reveal 
macroscopic  fracture  plane  inclinations.  Pepper's  Mine, 
located  near  the  Lajitas  station,  gives  access  to  the  subsur¬ 
face,  therefore  allowing  the  measure  of  fracture  plane  orien¬ 
tations  at  depth  with  the  use  of  a  Brunton  compass.  The 
infrared  satellite  photo  and  the  7  1/2  minute  Amarilla  Moun¬ 
tain  topographic  quadrangle  provides  the  means  to  construct  a 
lineation  map  representing  the  orientation,  size,  and  dis¬ 
tribution  of  macrofracture  traces  exposed  at  the  surface. 
Using  mylar  and  rapidograph  to  trace  visible  macrof racturc 
patterns  from  the  infrared  photo  produces  an  accurate  repre¬ 
sentation  of  fracture  sets.  Methods  discussed  in  Chapter  5 
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of  Compton's  1962  Manual  of  Field  Geology  suffice  in  deter¬ 
mining  scale  and  bearing  on  the  constructed  lineation  map. 
Upon  completion  of  the  lineation  map,  one  directly  measures 
macrofracture  separation  distances  and  orientations  using  a 
protractor  and  scale.  Because  fracture  orientation  varies 
with  geographical  location,  one  must  select  unit  volumes  rep¬ 
resentative  of  macrofracture  orientation  for  each  locality  of 
interest.  Rose  diagrams  facilitate  examination  of  dominant 
fracture  orientations  for  each  unit  volume  with  errors  in 
orientation  less  than  5  degrees. 

Determination  of  Microfracture  Distribution  and  Orienta¬ 
tion  .  Measuring  orientation  of  hairline  fractures  in  several 
oriented  in  situ  Santa  Elena  hand  specimens  yields  dominant 
microfracture  trends.  Using  a  Brunton  compass,  one  can  mea¬ 
sure  orientation  of  microfractures  within  3  degrees  before 
extracting  a  sample.  Rose  diagrams  document  the  spatial 
change  in  dominant  fracture  orientations. 

Collection  of  Oriented  Hand  Specimens.  From  each  sam¬ 
ple,  laboratory  work  requires  analysis  of  three  cylindrical 
plugs  measuring  2.54  cm  in  diameter  and  3  -  5  cm  in  length. 
Due  to  torsional  motion  in  the  coring  bit,  a  sample  riddled 
with  open  fractures  reduces  the  chance  of  recovering  plugs  of 
sufficient  length.  Hence,  one  must  choose  samples  large 
enough  and  with  as  few  fractures  as  possible.  Also,  plugs 
having  few  or  no  fractures  more  closely  reflect  the  nature  of 
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the  rock  matrix.  One  can  orient  samples  using  the  methods  of 
Prior  et  al .  (1987) . 

Observations 

There  is  a  correlation  between  microscopic  and  macro¬ 
scopic  fracture  orientation  and  distribution.  The  broken-up 
character  of  the  Santa  Elena  limestone  (Figure  2.3)  corre¬ 
sponds  to  dense  brittle  fracturing  observed  during  geologic 
reconnaissance  at  the  smaller  scale  in  outcrop. 

At  the  macroscopic  scale,  fracture  traces  unaffected  by 
topography  indicate  vertical  fracturing.  In  addition,  tops 
and  bottoms  of  massive  stratigraphic  layers  following  topo¬ 
graphic  contours  suggest  horizontal  bedding  (Billings,  1972). 
Three  distinct  fracture  systems  exist  (Figure  3.1).  However, 
rose  diagrams  from  areas  I,  Ila,  Ilb,  III,  and  IV  show  two 
general  macrofracture  trends  approximately  at  right  angles  to 
one  another--NNW  and  ENE  (Figure  3.2a-e).  Figure  3.2f  shows 
a  composite  of  fracture  orientations  for  all  five  areas. 
Fracture  separation  distance  ranges  from  35  -  40  meters. 

The  Santa  Elena's  "cobblestone"  (or  "blocky")  appearance 
suggests  subvertical  to  vertical  fracturing  and  horizontal 
parting  along  bedding  planes.  Also,  a  characteristic  cubic 
or  rectangular  hand  specimen  shape  implies  vertical  fractur¬ 
ing  and  horizontal  parting  at  the  microscopic  scale.  Two 
sets  of  fractures,  NNW  and  ENE,  dominate  at  the  microscopic 
scale.  Rose  diagrams  (Figure  3.3a-g)  represent  the  general 
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Figure  3.2a-e  Observed  strikes  of  macrof ractures  in 
unit  areas  a)  I,  b)  Ila,  c)  lib,  d)  III,  and  e)  IV. 
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Macrofracture  Orientations 
Units  I,  Ila,  Ilb,  III,  and  IV 


=  200 


Figure  3.2f  Cummulative  representation  of  observed 
macrofracture  strikes  in  areas  I  through  IV. 


Figure  3 . 3a-g  Observed  strikes  of  microfrac¬ 
tures  in  samples  a)  VLS-1,  b)  VLS-2,  c)  VLS-3, 
d)  VLS-5,  e)  VLS-6,  f)  VLS-7,  and  g)  VLS-8 . 
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trends  of  hairline  fracture  planes  in  each  sample.  Figure 
3 . 3h  shows  a  composite  of  all  fracture  plane  strikes  at  the 
microscopic  scale.  Thick  caliche  cover  on  sample  VLS-4  and 
VLS-8  camouflaged  a  dense  distribution  of  hairline  fractures 
visible  only  on  cutting  the  samples.  Consequently,  only  a 
few  fractures  could  be  measured  from  sample  VLS-8.  No 
fractures  were  measured  from  sample  VLS-4.  The  fracture  sep¬ 
aration  distance  at  the  microscopic  scale  ranges  from  2.5  - 
30.5  cm . 

Pepper's  Mine,  once  mined  for  its  mercury,  provides 
access  to  the  subsurface  (Figure  3.4) .  A  dense  distribution 
of  vertical  fractures  lines  the  tunnel  walls.  The  tunnels 
follow  fractures  oriented  with  either  of  the  two  general 
trends  measured  from  the  lineation  map  and  oriented  hand 
specimens  (Figure  3.5)  .  Excluding  frequent  horizontal  part¬ 
ing  along  bedding  planes,  most  fractures  have  subvertical  to 
vertical  inclinations.  Although  most  bedding  is  horizontal, 
in  a  few  areas,  undulating  laminated  sandy  limestone  beds 
ranging  from  10  -  15  cm  in  thickness  and  massive  "pinched 
out"  beds  cause  a  variation  in  bedding  dip  from  14°  to  31°. 
Fractures  and  partings  between  bedding  planes  are  healed  and 
partially  healed  with  sparry  calcite  and  cinnibar  fillings. 
Fillings  separating  beds  range  from  7  -  30  cm  in  thickness, 
whereas  fracture  plane  fillings  range  from  0  -  7  cm  in  thick¬ 
ness.  Because  of  possible  blast  induced  effects,  fracture 
separation  distance  in  the  mine  was  not  estimated. 
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Microfracture  Orientations 
Samples  VLS-1,  VLS-2,  VLS-3,  VLS-5,  VLS-6 
VLS-7  and  VLS-8 


Figure  3 . 3h  Cummulative  representation  of  observed 
microfracture  strikes  in  samples  VLS-1,  VLS-2,  VLS-3, 
VLS-5,  VLS-6,  VLS-7  and  VLS-8. 
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SKETCH  OF  PEPPER'S  MINE 
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Fracture  Orientations 
Mine 


The  Santa  Elena  is  comprised  primarily  of  calcite  and 
secondary  clay,  silica,  and  trace  amounts  of  pyrite.  Samples 
range  in  color  from  buff-  to  yellow-  to  violet-white  and  in 
some  cases  two  to  three  colors  intermingle,  giving  a  marbled 
appearance.  Bioturbation  caused  this  intermingling.  The 
texture  is  soft  (H  =  3  -  3  1/2),  fine-grained  and  massive — a 
few  samples  look  chalky.  Fossil  content  varies  from  approxi¬ 
mately  0  to  5  percent.  Observed  Nerinied  (class  gastropoda), 
Caprinid  (class  pelecypoda) ,  ammonite  and  foraminifera  fos¬ 
sils  show  sparry  calcite  replacement.  However,  silicifica- 
tion  was  observed  in  several  rudistids  (phylum  echinoder- 
mata) .  Fresh  cut  surfaces  revealing  holes  less  than  0 . 5  mm 
coupled  with  quick  evaporation  in  some  water  saturated 
samples,  suggests  a  connected  porosity  range  of  0  to  10 
percent.  Healed  and  partially  healed  fractures  present  in 
hand  specimens  have  sparry  calcite  fillings. 

In  summary,  three  fracture  sets  exist  in  the  Santa  Elena 
lirr^stone  at  the  Lajitas,  Texas  site:  vertical  NNW  and  ENE 
tre.nai.ng  orthogonal  sets,  and  horizontal  parting  along  bed¬ 
ding  planes.  This  fracturing  occurred  both  on  a  microscopic 
a.nd  macroscopic  scale.  Based  on  field  observations,  the 
Santa  Elena  limestone  is  composed  primarily  of  calcite  and 
varies  in  porosity  from  sample  to  sample. 
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IV.  Laboratory  Measurements 


Determining  microscopic  fracture  orientation  and  seismic 
velocity  in  hand  specimens  in  the  laboratory  aids  in  measur¬ 
ing  microscopic  inhomogeneities  (such  as  vuggy  porosity,  high 
clay  content,  fractures,  etc.)  that  may  significantly  affect 
propagating  waves  through  material  in  situ.  Testing  the 
Santa  Elena's  degree  of  homogeneity  in  the  laboratory  entails 
methods  similar  to  those  used  by  Sorrells  (1961)  .  The 
methodology  and  laboratory  results  for  measurements  of  P-  and 
S-wave  velocity,  dry  and  saturated  bul)<  roc)c  density,  grain 
density,  porosity,  microscopic  fracture  and  pore  geometry, 
and  mineralogical  composition  follow  in  the  next  section. 

Note  that  the  total  porosity  <[)  of  a  bullc  medium  equals  the 
sum  of  its  total  crack  porosity  <{)c  and  total  pore  porosity  Op; 
where  Oc  is  the  total  void  space  between  open  fractures  arid  Op 
is  the  total  spheroidal  to  ellipsoidal  void  space  within  a 
given  unit  volume . 


Methodology 

Porosity/Density  Measurement.  Spatial  changes  in  poros¬ 
ity  and  composition  can  cause  variations  in  bulk  rock  veloci¬ 
ties  of  as  m'uch  as  20  percent.  Heterogeneities  sucii  as  these 
can  effect  velocities  measured  in  situ.  Without  kr'O.w.fviae  of 
bulk  rock  porosity  and  compositional  variations,  one  :tay 
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falsely  conclude  that  azimuthal  anisotropy  exists. 

Laboratory  porosity  and  dry  bulk  rock  and  grain  density 
measurements  aid  in  determining  variation  in  composition  and 
porosity  between  samples. 

Determination  of  the  mineral  grain  density,  saturated 
and  dry  bulk  rock  densities  and  porosity  requires  the 
measurement  of  sample  weight  parameters .  These  parameters 
include  dry  weight,  w^,  suspended  weight,  W5J,  and  saturated 
weight,  Wsa  .  The  procedure  for  these  measurements  utilizes 
the  techniques  developed  by  Archimedes  in  210  B.C.  (Tipler, 
1982) .  To  obtain  dry  weight,  one  must  first  extract  all 
moisture  from  plug  samples  using  a  vacuum  oven  set  at 
approximately  60°C.  Leaving  samples  in  the  vacuum  oven  for 
24  hours  and  then  allowing  the  samples  to  cool  20  minutes 
before  weighing  on  the  Mettler  PE  3600  yields  dry  weight 
measurements  with  errors  less  than  0.01  grams.  Prior  to 
measuring  saturated  and  suspended  weights,  one  must  saturate 
the  samples  using  a  pressure  vessel  (Figure  4.1)  . 

Obtaining  saturated  and  suspended  weights  entails  weigh 
ing  samples  both  resting  on  the  bottom  of  and  suspended  in  a 
beaker  containing  "saturation  fluid";  in  this  case,  water. 
The  beaker  rested  on  a  Mettler  PE  3600  balance  capable  of  re 
solving  weights  to  ±0.01  grams.  To  determine  the  saturated 
weight,  wg^ : 

(4.1)  Wsa  =  “  '''.-■-3 
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Figure  4.1  Pressure  vessel  used  to  saturate  plug 
samples  with  water. 


where  Ws+f^b  equals  the  weight  of  the  fluid,  beaker,  and  a 
saturated  sample  resting  fully  immersed  on  the  bottom  of  the 
beaker;  and  wf*b  equals  the  weight  of  the  beaker  and  fluid. 

To  determine  the  suspended  weight,  wg^ : 

(^•2)  Wso  =  Wsus  +  F*B  “  Wr  +  B 

where  wsus+f  +  b  equals  the  weight  of  the  fluid,  beaker,  and  a 
fully  immersed  saturated  sample  suspended  from  a  string. 
Measuring  wp-^g  before  lowering  each  sample  into  the  fluid- 
filled  beaker  is  important.  Due  to  surface  tension,  the  vol 
ume  of  fluid  in  the  beaker  decreases  with  each  sample 
removal,  thus  wp^g  changes  with  each  progressive  set  of  satu¬ 
rated  and  suspended  weight  measurements.  Also,  using  a  fine 
string  such  as  fishing  line  to  slowly  lower  sa.mples  into  the 
beaker  eliminates  fluid  loss  due  to  splashing. 

Measurements  of  the  dry  weight  wg,  wg-j,  and  wsa  are  suff 
cient  to  calculate  connected  porosity  and  density  for  each 
sample.  Taking  the  ratio  of  pore  volume  to  total  volume, 
Vp/Vx,  yields  fractional  porosity,  0,  where: 

(4.3)  Vp  =  (wsA  -  wg)  /p. 

and 

(4.4)  Vj  =  Wg-/ p,  . 

Thus,  with  saturation  fluid  density,  P;,  equivalent  to  1.00 
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gm/cc  for  water: 


(4.5)  *t*  “  Vp/V^  —  (WsA  “  wq)/wsu. 

For  dry  and  saturated  rock  densities,  Pdry  and  Psac : 

(4.6)  Pdry  ~  Wq/Vx  =  (Wq/ Wsu)  / Pj  , 
and 

(4.7)  Psat  ~  '^Sa/^T  ~  (  Wsa/Wsu  ) /p  J  . 

The  grain  density,  Pgrain>  equals: 

(4.8)  pgrain  =  Wq/Vs  =  Wq/  [Pf  ( Wsu  +  Wr.-WsA )  ] 

where  Vs  represents  the  solid  volume  of  a  sample,  the 
difference  between  the  total  plug  volume  and  pore  volume: 

(4.9)  Vx  -  Vp  =  (Wsu  +  Wq-Wsa)  /P;  . 

This  method  gives  accurate  estimates  (±0.005)  of  dry 
bulk  and  saturated  rock  densities,  grain  density,  and  poros¬ 
ity  provided  unconnected  void  space  does  not  exist  in  concen¬ 
trations  greater  than  0.5  percent.  For  instance,  a  dry  sam¬ 
ple  of  pure  calcite  equal  to  a  bulk  volume  of  19.306  cc  with 
void  space  equal  in  concentration  to  0.06  for  air-filled  con¬ 
nected  porosity  and  0.03  for  isolated  water-filled  porosity 
weighs  48.3654  g  and  equals  2.5052  g/cc  in  density: 
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(4.10) 


Pclry  =  (l-<l>)pS-^^‘t)„)p„, 

=  (l  -  0.09)  ^2.12  +  (0.03)  (l.OO 

=  2.5052 


cc 


cc 


and 


(4.ii; 


=  (Pd) 


V 


T  ! 


=  [2.5052  19.306  cc  =  48.3654  g  , 

where  the  total  porosity  and  isolated  porosity  is  (})  and  0^ 
respectively  and  the  density  of  water  is  pw .  The  saturated 
weight  equals 


(4.12) 


WSA  =  (U-  0)  Ps+(0pPw))  ^ 
=  49.5238  g  . 


Applying  the  rules  of  Archimedes: 


(4.13) 


(4.14: 


P  qrain 


Wsu  =  Vt  P„=  19 . 306  cc  , 


48.3654  q 


1.00  g/cc  (19.306  g  +  49.3654  g  -  49.5238  g) 
=  2.5259  g/cc 


and 


(4.15) 


0  =  49.5238  g  -  48.3654  g  ^  q  q. 
19.306  g 
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Thus,  the  methods  of  Archimedes  do  not  account  for  isolated 
pore  porosity  and  consequently  underestimate  grain  density 
when  isolated  porosity  exists. 

Table  4.1  lists  density  and  porosity  measurements  for 
eacn  of  three  plugs  extracced  f-Ciii  si::  of  oho  alght  orient^^d 
Santa  Elena  hand  specimens.  (Plugs  of  sufficient  length  could 
not  be  obtained  from  densely  fractured  samples  VLS-4  and  VLS- 
8.)  Bulk  dry  rock  densities  range  from  2.337  gm/cc  to  2.611 
gm/cc,  whereas  grain  densities  range  from  2.667  gm/cc  to 
2.721  gm/cc.  Calculations  omitting  the  two  grain  density 
extremes  yield  a  mean  grain  density  of  2.70  gm/cc  plus  or 
minus  a  0.03  maximum  deviation.  Thus,  with  grain  densities 
measuring  no  less  than  98  percent  of  pure  calcite’s  2.72 
gm/cc  density,  differences  observed  between  dry  bulk  and 
grain  density  in  Table  4.1  are  most  likely  due  to  connected 
pore  space  or  air-filled  microfractures  and  small  concentra¬ 
tions  (less  than  3  percent)  of  isolated  pore  space.  Grain 
densities  less  than  2.70  gm/cc  imply  the  presence  of  sec¬ 
ondary  minerals  and/or  unaccounted  isolated  pore  space.  Per¬ 
cent  porosity  correlates  well  with  changes  between  dry  bulk 
and  grain  densities  from  plug  to  plug. 

Determination  of  Mineraloaical  Composition.  Pure  lime¬ 
stones  with  isolated  pore  space  yield  grain  densities  less 
than  that  of  pure  calcite  (Pcalcite  =  2.72  g/cc) .  Grain  den¬ 
sities  less  than  2.72  g/cc  (by  approximately  2  percent  or 
more)  also  result  when  limestones  contain  secondary  minerals 
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in  concentrations  approximately  greater  than  10  percent. 
Knowledge  of  mineralogical  composition  therefore  facilitates 
estimates  of  total  porosity  in  hand  specimens. 

X-ray  diffraction  analysis  effectively  identifies  miner¬ 
als  and  their  corresponding  relative  abundances  present  in 
rock  specimens.  This  type  of  analysis  examines  diffraction 
effects  of  x-r^,vs  incident  on  crystal  lattice  structures. 

Determining  mine ral : g ica 1  composition  using  x-ray 
diffraction  techniques  requires  four  basic  steps:  1)  sample 
preparation  and  estimation  of  coarse  to  fine  particle  weight 
percents,  2)  measurement  of  diffracted  angles,  3)  application 
of  Bragg's  law,  end  4)  identification  of  minerals  and  estima¬ 
tion  of  their  relative  abundances. 

To  avoid  grinding  distinct  grains  during  sample  prepara¬ 
tion,  1  -  2  mm  disks  sliced  from  each  plug  sample  were 
crushed  using  a  mortar  and  pestle.  Determining  the  relative 
weight  percents  of  sand  and  silt  and  clay  size  particles 
necessitated  the  measure  of  the  total  crushed  sample  weight, 
wtot,  and  the  weight  of  sand  and  silt,  Wss^  for  each  sample. 
Measurement  of  Wgs  required  the  separation  of  sand  and  silt 
sized  particles  from  clay  size  particles  less  than  4  microns. 
After  soaking  each  crushed  sample  in  de-ionized  water  and 
sodium  pyrophosphate  to  minimize  flocculation  problems,  a 
sonic  probe  aided  in  removing  clay  size  particles  (or  grain 
coatings)  irom  host  grains.  Centrifuging  and  decanting  the 
liquid  separated  sand  and  jilt  grains  from  clay  particles  in 
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suspension.  Drying  and  weighing  residues  left  in  the  bottoms 
of  the  centrifuge  tubes  yielded  weight  percents  of  coars^^  and 
fine  particles: 

(4.16)  100  (wss/w-tot)  “  %  sand  and  silt  size  particles 

(4.17)  100  -  [  100  (wss/wtot)  ]  =  %  clay  size  particles 

A  micronizing  machine  ground  the  sand  and  silt  residues  to 
white  flour;  randomly  orienting  grains  and  consequently  pro¬ 
ducing  good  diffraction  characteristics.  Measuring  the 
change  in  d  spacing  from  a  dry  to  saturated  state  helped  in 
identifying  the  clay  type.  To  measure  d  spacings  along  the  c 
axis  in  clays,  clay  particles  left  in  suspension  were  fun¬ 
nelled  through  filters.  (Due  to  their  platy  nature,  the 
particles  tended  to  settle  on  flat  faces  )  The  d  spacing  in 
clays  equals  the  distance  between  indefinite  '“xtended  s^^eets 
of  Si04  tetrahedra;  where  three  of  the  four  ox^  ens  in  each 
Si04  tetrahedron  are  shared  with  neighboring  tetrahedra, 
leading  to  a  ratio  of  Si:0  =  2:5.  Each  sheet,  if 
undistorted,  has  a  hexagonal  symmetry  with  the  c  axis 
perpendicular  to  the  sheets  of  tetrahedra  (Hurlburt  and 
Klein,  1977)  .  Because  some  clays  swell  or  expand  when 
saturated  with  fluid,  d  spacings  can  become  enlarged.  The 
degree  of  d  spacing  enlargement  depends  on  the  type  of  clay. 
Hence,  diffraction  patterns  were  obtained  for  oriented  clay 
specimens  both  in  an  air  and  liquid  saturated  state. 
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Ethylene  glycol,  a  fluid  that  expands  and  stabilizes  d 
spacing,  was  used  to  saturate  clays  (Klug  and  Alexander, 

1974)  . 

X-ray  diffraction  analysis  revealed  the  presence  of  sec¬ 
ondary  minerals  in  abundance  of  2.5  to  8.8  percent  (Table 
4.2) .  These  secondary  minerals  include  quartz,  ranging  in 
quantity  from  0.8  to  6.1  percent,  and  clays,  ranging  from  1.3 
to  2.7  percent.  Relative  percentages  for  quartz,  calcite, 
and  total  clay  minerals  have  a  2  percent  error  factor; 
whereas  the  percentage  error  present  in  relative  clay 
abundances,  such  as  smectite,  illite,  and  )caolinite, 
approaches  20  percent.  Calcite,  ranging  from  91.2  to  97.5 
percent  (in  relative  abundance  from  plug  to  plug),  is 
obviously  the  primary  constituent.  Therefore,  results  from 
the  X-ray  diffraction  analysis  support  the  inferences  made 
from,  the  porosity/density  measurements.  Differences  observed 
between  dry  bul)c  and  grain  density  result  primarily  from  open 
pore  space  or  microfractures,  rather  than  substantial 
secondary  mineral  concentrations. 

Velocity  Measurement.  Detecting  and  quantifying  veloc¬ 
ity  anisotropy  requires  knowledge  of  either  P-  or  S-wave 
velocity.  For  P-waves,  defining  anisotropy  resulting  from 
wave  propagation  through  a  fractured  solid  requires  measure¬ 
ment  of  compressional  velocity  along  two  or  more  physically 
different  propagation  paths.  Factors  other  than  anisotropy, 
such  as  variation  in  rock  composition,  pore  fluid  saturation, 
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TABLE 
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otal  100.0%  IM.0%  100.0%  100.0%  100.0%  100.0%  100.0%  100.0% 


etc.,  also  influence  P-wave  velocity  more  than  S-wave 
velocity.  Shear  waves  give  more  information  about  aligned 
cracks  than  do  P-waves.  Shear  waves  in  anisotropic  media 
split  into  two  orthogonal  polarized  components  which  travel 
along  identical  propagation  paths.  These  two  components 
propagate  at  different  velocities  and  separate  in  tim.e.  In  a 
fractured  solid,  the  shear  wave  component  polarized  parallel 
to  the  fracture  plane  strikes  arrive  first  at  the  receiving 
end  (Crampin,  1985)  .  Thus  measuring  both  P-  and  S-wave 
velocity,  with  shear-wave  displacements  polarized  both 
parallel  and  perpendicular  to  observed  fracture  trends  in 
hand  specimens  and  outcrop,  aids  in  defining  pore  geometry 
(spheroids  versus  fractures) .  Also,  Vp/Vs  gives  information 
conoerning  mineralogical  composition,  in  addition  to  that 
obtained  from  porosity/density  measurements  and  XRD  analysis 
(Wilkens  et  al.,  1984). 

Velocity  was  measured  for  3  -  5  cm  cylindrical  plugs 
2.54  cm  in  diameter;  as  plugs  longer  than  5  cm  vibrate  and 
generate  tube  waves  that  produce  spurious  results.  Error  in 
travel  time  measurements  increases  from  ±0.03  km/s  to  ±0.06 
km/s  when  plugs  are  less  than  2.54  cm  in  length. 

Ends  were  trimmed  and  ground  flat  parallel  to  one 
another  within  0.0015  cm  for  three  plugs  extracted  from  each 
of  the  six  oriented  Santa  Elena  rock  samples .  The  orienta¬ 
tion  of  each  plug  axis  parallels  poles  to  fractures  or  bed¬ 
ding.  The  axes  of  plugs  Y  and  Z  parallel  northwest  and  north- 
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east  fracture  plane  normals  respectively,  whereas  the  X  axis 
plug  is  vertical  (Figure  4.2) .  Orientation  of  each  plug  axis 
comes  within  10  degrees  of  fracture  and  bedding  plane  nor¬ 
mals. 

In  the  laboratory  the  velocity  V  is  determined  from  the 
sample  length  L  and  the  transit  time  AT; 

(4.18)  V  =  L/AT 

An  ultrasonic  pulse  technique  developed  by  Simmons  (1965) 
measures  travel  times  with  errors  less  than  0.05 
microseconds . 

The  travel  time  of  an  elastic  wave  through  a  set  of 
transmitting  and  receiving  transducers  equals  the  time  delay 
between  a  pulse  input  and  signal  output  (minus  any  correc¬ 
tions)  .  Because  waves  must  pass  through  a  coupling  resin  and 
a  transducer  facing  material,  designed  to  reduce  ringing, 
static  time  corrections  must  be  subtracted  from  raw  measured 
transit  times.  Correction  times  equal  0.65  and  1.00  |ls  for 

P-  and  S-waves  respectively. 

With  all  measurements  made  under  ambient  conditions 
(i.e.,  bench-top)  with  no  saturating  liquids  or  confining 
pressure,  t.he  application  of  a  uniaxial  load  of  appro.ximate iy 
5-10  bars,  with  the  use  of  parallel  clamps,  improved  cou¬ 
pling  between  plug  end  and  transducer  surfaces,  and  thus 
yielded  good  P-  and  S-wave  data. 
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Figure  4.2  Illustration  showing  plug  axes 
orientations . 
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Two  switches  in  the  transducer  electrical  circuit 


allowed  for  the  change  from  longitudinal  displacement  to 
either  of  two  shear  wave  polarizations.  The  propagation 
paths  and  particle  displacement  orientations  of  these  three 
waves  differ  slightly,  as  shown  in  Figures  4.3  and  4.4. 
Transmitting  one  wave  at  a  time  substantially  reduced  P-  to 
S-  and  S-  to  P-wave  conversions  due  to  reflections  off  the 
plug  wall . 

The  frequency  band  width  ranged  from  1.5  -  2.75  MHz  with 
wavelengths  averaging  approximately  1-2  millimeters.  Each 
transducer  disk  used  in  this  study  had  a  free  resonant 
frequency  of  1  MHz. 

The  velocity  data  support  the  inferences  made  from  the 
porosity/density  measurements  and  XRD  analysis  (Table  4.3). 
Shear  wave  velocities  vary  azimuthally  less  than  2.00  percent 
(and  in  most  cases  less  than  0.35  percent)  within  each 
individual  plug.  This  2.00  percent  variation  is  well  within 
the  bounds  of  experimental  error,  hence  azimuthal  anisotropy 
either  does  not  exist  or  can  not  be  resolved  (solely 
comparing  fractional  porosity  to  wave  velocity)  in  the 
laboratory  samples.  However,  shear  velocities  vary  sub¬ 
stantially  from  plug  to  plug — in  the  most  extreme  case  as 
much  as  12  percent.  Variatloii  of  compressional  wave  veloci¬ 
ties  from  plug  to  plug  ranges  from  0.7  -  15.8  percent.  The 
higher  values  of  measured  connected  porosity  correlate  with 
the  lower  P-  and  S-wave  velocity  measurements.  Appendix  B 
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Figure  4.3  Configuration  of  upper  and  lower 
transducer  sets. 
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compressional  and  shear  waves  generated  by  transducers 
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LABORATORY  VELOCITY  AND  POROSITY  MEASUREMENTS 
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gives  a  representative  sampling  of  the  waveforms  from  which 
travel  times  were  selected. 

Figures  4.5  and  4.6  illustrate  the  changes  in  P-  and  S- 
wave  velocity  with  varying  porosity.  Scatter  exceeds  the 
limits  of  experimental  error  for  porosities  greater  than  5 
percent.  Distribution  and  variation  of  pore  geometries  and 
changes  in  clay  and  silica  content  less  than  a  few  percent: 
may  contribute  to  this  scatter. 

The  P-  and  S-wave  matrix  velocities,  obtained  from  the 
zero  intercepts,  equal  6.04  km/s  and  3.23  km/s  respectively. 
These  matrix  velocities  fall  below  Dandekar's  (1968)  matrix 
Vp  and  Vs  measurements  of  6.53  km/s  and  3.36  km/s  for  pure 
calcite.  Figure  4.7  shows  the  Vp/Vs  versus  porosity  for  each 
plug  cf  Ccnta  Elena  limestone.  Again  the  zero  intercept,  of 
the  linear  curve  best  fit  to  the  velocity  data,  falls 
slightly  below  Dandekar's  (1968)  Vp/Vs  measure  of  1.94.  In 
contrast  to  calcite,  the  matrix  velocities  and  Vp/Vs  pure 

quartz  equal  6.05  km/s  (Vp) ,  4.09  km/s  (Vs),  and  1.48  (Vp/Vs) 
(McSkimin  et  al . ,  1965) . 

The  influence  of  composition,  porosity,  and  pore  geom.e- 
try  must  now  be  quantified  in  order  to  interpret  the  labora¬ 
tory  seismic  data.  For  a  set  of  siliceous  limestones, 

Wilkens  et  al.  (1984)  measured  Vp  and  Vs/  determined  bulk  den¬ 
sity,  effective  porosity,  and  carbonate  content;  a.nd  observed 
pore  and  fracture  distributions  using  a  scanning  electron 
microscope.  They  showed  that  composition  is  more 
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P-WAVE  VELOCrTY  VS.  POROSITY 
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Vp/Vs  Ratio  VS.  Porosity 


POROSITY  (») 

Figure  4.7  Vp/Vs  ratio  versus  porosity  for  each 
plug  of  Santa  Elena  limestone. 
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important  than  porosity  and  pore  geometry  for  Vp/Vs-  Accord¬ 
ing  to  Wilkens  et  al.  (1984),  the  effect  of  composition  is 
best  ascertained  from  data  obtained  at  high  pressures  where 
microfractures  are  closed. 

At  1  kbar,  both  Vp  and  Vs  increase  with  increasing  bulk 
density  and  decreasing  porosity  in  composit ionally  homoge¬ 
neous  sample  sets.  But,  for  siliceous  limestones,  Vp 
increases  monotonically  and  Vs  remains  invariant  because 
quartz,  although  less  dense  than  calcite,  has  a  higher  shear 
wave  velocity.  (Vp/Vs  for  quartz  equals  1.5;  whereas  for 
common  phase  minerals  in  sedimentary  rocks  such  as  calcite, 
dolomite,  and  feldspar,  Vp/Vs  equals  or  falls  between  1.8  and 

2.0.)  Increasing  pressure  (closing  microfractures)  has  lit¬ 
tle  effect  on  Vp/Vs  between  0.01  and  1.0  kbar.  Changing 
porosity  by  ±  5  percent  in  silica-rich  limestones  decreases 
Vp/Vs  by  approximately  5  percent,  but  as  quartz  content 
decreases  and  calcite  increases,  Vp/Vs  increases  by  20  per¬ 
cent.  That  is,  total  porosity  and  pore  geometry  cause  varia¬ 
tions  of  approximately  0.1  in  Vp/Vs,  whereas  compositional 
variation  equals  0.4.  Because  the  Santa  Elena's  Vp/Vs  of  1.87 
comes  so  close  to  the  calcite  ratio,  the  XRD  measurements  of 
1.3  to  2.6  percent  in  clay  content  suggest  that  clay  probably 
plays  a  role  in  lowering  the  Santa  Elena's  P-  and  S-wave 
matrix  velocities,  as  concentrations  of  only  3  percent  can 
lower  longitudinal  wave  velocities  (Sorrells,  1961). 
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For  pure  limestones,  pore  aspect  ratios  determined  with 
the  theory  of  Toksoz  et  al.  (1976)  (for  seismic  velocities) 
agree  with  the  Wilkens  et  al .  (1984)  electron  microscope 

observations.  Based  on  this  earlier  success,  laboratory  mea¬ 
surements  of  the  Santa  Elena's  P-  and  S-wave  velocity  normal¬ 
ized  to  the  P-  and  S-wave  matrix  velocities  of  6.05  km/ s  and 
3.23  km/ s  respectively,  are  plotted,  on  the  theoretical 
curves  of  Toksoz  et  al.  (1976)  for  equations  (4.22)  and 
(4.23),  as  a  function  of  aspect  ratio  and  gas-saturated 
porosity  in  Figures  4.8  and  4.9.  Aspect  ratio,  defined  as 
the  ratio  of  the  minor  to  major  semiaxis  of  a  spheroid,  esti¬ 
mates  pore  geometry.  Data  points  for  P-waves  range  in  aspect 
ratio  from  approximately  0.1  to  1.0.  With  the  exclusion  of  6 
out  of  32  data  points  greater  than  1.0,  S-wave  data  points 
also  range  from  approximately  0.1  to  1.0.  For  all  S-wave 
data  points  to  plot  at  or  berow  1.0,  the  normalizing  velocity 
must  equal  3.30  km/s;  hence  the  S-wave  intercept  velocity 
chosen  from  Figure  4.6  may  be  low  by  0.07  km/s.  Experimental 
error  may  be  another  contributing  factor  to  the  high  aspect 
ratios.  Figures  4.8  and  4.9  indicate  that  pore  spaces, 
ellipsoidal  to  spherical  in  shape,  primarily  contribute  to 
the  scatter  observed  in  Figures  4.6  and  4.7. 

In  order  to  utilize  results  such  as  those  in  Figure  4.8 
and  4.9,  it  is  necessary  to  investigate  the  foundation  of  the 
Toksoz  et  al.  (1976)  model.  Toksoz  et  al.  (1976)  demonstrate 
the  effects  of  inclusion  shapes  on  the  velocities  of  a 
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NORMALIZED  S-WAVE  VELOCITY  VS.  POROSITY 
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POROS  lY  (») 

Figure  4.9  Theoretical  shear  wave  velocity-porosity  relationships  for  differing 
aspect  ratios  versus  data  from  this  study.  The  ratio  of  the  minor  to  major 
semiaxis  of  a  spheroid  defines  the  aspect  ratio  a. 


composite  medium  by  comparing  velocities  and  reflection  coef¬ 
ficients  predicted  by  an  idealized  model  with  published  labo¬ 
ratory  results.  Their  model  defines  a  porous  rock  in  terms 
of  a  solid  elastic  matrix,  randomly  distributed  pores,  and 
saturating  fluids  (gas,  oil,  or  water) .  Assuming  wavelengths 
long  in  comparison  to  pore  size,  and  utilizing  the  formulas 
of  Kuster  and  Toksoz  (1974)  based  on  scattering  theory 
enables  the  approximation  of  a  whole  rock  in  terms  of  an 
equivalent  homogeneous  medium  with  some  effective  elastic 
coefficients.  Toksoz  et  al.  (1976)  use  these  equations  for 
calculating  velocities  and  extend  the  theory  to  cover  mixed- 
fluid  (i.e.,  gas-water,  oil-water)  saturation  for  a  spectrum 
of  shapes.  They  also  derive  expressions  relating  pressure  to 
seismic  velocities  in  a  porous  and  saturated  medium. 

Specifying  the  medium  properties  in  terms  of  the  bulk 
modulus  K,  shear  modulus  |i,  and  density  p 

K-K  1  K'-K  ^ 

~  o  •  ^  •  ‘  a  JJ  ' 

/A  T  q\  3K+4p.  3K+4)J. 


(4.20) 


p-p 

T  — .  ^  T 

6p  (k+2p)+  p(9K+8p 

-c 

j  25p(3K+4p) 

and 


(4.211 


p  =  p  (1-c)  +  cp  '  , 
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where  c  is  the  volume  concentration  of  the  inclusions,  I'iijj 
and  Tijxj  are  scalar  quantities,  functions  of  K,  4,  K',  4', 
and  the  aspect  ratio,  a,  of  the  inclusions.  The  ratio  of  the 
minor  to  major  semiaxis  of  a  spheroid  defines  th^  aspect 
ratio.  A  fracture’s  aspect  ratio  is  a  measure  of  its  flat¬ 
ness.  Unprimed  quantities  (K,  4,  p)  and  primed  quantities 
(K',  4',  p')  refer  to  the  matrix  and  inclusions  respectively 

whereas  "tilded"  quantities  equal  effective  properties  of  the 
composite  medium.  (See  Appendix  C  for  expressions  of  T^ajj 
and  Txjsj.) 

Using  a  rock  model  derived  on  the  basis  of  laboratory 
data  in  the  equations  for  calculating  velocities,  Toksdz  et 
al.  (1976)  reach  several  conclusions.  The  conclusions  most 
important  to  this  study  are  1)  for  a  given  matrix  the  compos¬ 
ite  elastic  moduli  and  seismic  velocities  of  rocks  decrease 
with  increasing  porosity,  2)  compressional  velocities  are 
affected  more  by  properties  of  saturating  fluids  than  shear 
velocities,  and  3)  for  a  given  pore  concentration,  flatter 
(thinner)  pores  affect  velocities  substantially  more  than 
rounder  or  spherical  pores.  In  fact,  the  presence  of  pores 
with  aspect  ratios  less  than  0.01  in  concentrations  of  less 
than  a  percent  can  decrease  velocities  by  as  much  as  20 
percent . 

Typically  a  whole  rock  contains  pore  shapes  ranging  from 
nearly  equidimens iona 1  to  very  flat  thin  spaces.  Theoreti¬ 
cally,  this  variation  can  be  represented  in  terms  of  a  spec- 
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trum  of  aspect  ratios;  where  spheres  and  rounded  spheroids 
approximate  equidimensional  and  vugular  pores,  and  ellipsoids 
of  low  aspect  ratio  represent  fractures  or  flat  pores.  Two 
or  more  crack  populations  comprise  a  spectrum  of,  aspect 
ratios;  that  is,  cracks  in  population  "A"  have  aspect  ratios 
equal  to  "a"  and  those  in  population  "B"  have  aspect  ratios 
equal  to  "b" .  As  air  saturated  pores  of  small  aspect  ratio 
cause  a  greater  change  in  compressional  wave  velocities  than 
shear  velocities,  determining  the  effects  of  total  porosity 
on  seismic  velocities  in  a  whole  rock  requires  one  to  examine 
the  effects  of  each  population  individually.  Generalizing 
the  theoretical  formulas  given  by  equations  (4.19)  to  (4.21) 
for  the  case  of  mixed  aspect  ratios,  the  effective  bulk  and 
shear  modulus  equal 


(4.22) 
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3K+4^ 


1 

3  3K+4|1 


and 


(4.23) 


-  la  iaj 

6|a.(K+2|l)+  (l(9K+8|l) 

2  5^i(3K+4^i) 

where  c  {(X^)  equals  the  concentration  of  pores  with  aspect 
ratio  ttm  and  the  total  porosity  in  this  case  is 
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(4.24) 


4)= 

m=l 

Assuming  "noninteraction"  between  spheroidal  pores  and  frac¬ 
tures  (as  this  theory  does)  imposes  the  restriction 
[c  (OCm)  ] /Otm  <  1  in  the  application  of  equations  (4.22)  to 

(4.24)  (Toksciz  et  al  .  ,  1976). 

As  the  seismic  data  in  Figure  4.8  and  4.9  indicate,  the 
application  of  the  Toksdz  model  predicts  aspect  ratios 
between  0.1  to  1.0.  This  geometry  prediction  is  next  tested 
by  crack  section  examination. 

Examination  of  Microstructure.  Pore  geometry  can  be 
examined  in  crack  sections  under  the  microscope.  Examination 
of  crack  sections  60  -  150  microns  in  thickness  enables  a 
more  detailed  study  of  microstructure  than  standard  petro¬ 
graphic  thin  section  analysis.  Two  to  five  times  the  thick¬ 
ness  of  thin  sections,  crack  sections  produce  fewer  fractures 
during  the  sectioning  process  and  provide  three  dimensional 
views  of  microstructure.  Under  a  binocular  microscope,  the 
greater  thickness  of  section  creates  several  focal  planes, 
thus  allowing  the  tracking  of  microfractures  and  pore  spaces 
through  the  section.  Crack  sections  allow  the  recognition  of 
structural  features  not  detectable  in  thin  section  (Wilkens 
et  al .  ,  1984). 

Crack  sections  cut  to  a  few  microns  greater  than  the 
desired  thickness  on  a  thin  section  cut-off  saw  began  the 
sectioning  process.  Impregnation  of  sections  with  a  low  vis- 
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cosity  blue-stained  epoxy  in  a  vacuum  chamber  defined  frac¬ 
ture  and  pore  geometry  more  sharply.  Excess  epoxy  was  planed 
off  after  drying  and  sections  were  polished  smooth  on  a  thin 
section  grinder. 

Crack  section  analysis  of  the  Santa  Elena  limestone  sup¬ 
ports  the  inferences  and  conclusions  made  from  the  velocity, 
XRD,  and  porosity/density  data.  Two  of  five  sections  cut 
from  plug  ends,  XI  and  X3,  contain  several  fractures  healed 
with  sparry  calcite.  Sections  from  X5  and  Z7  contain  no 
fractures  and  only  one  section,  Y7,  contains  an  open  frac¬ 
ture.  Calcite  recrystallization  took  place,  especially  in 
sections  filled  with  masses  of  crushed  fossil  debris.  Pores 
are  vugular  and  commonly  isolated  from  one  another.  Isolated 
pore  space  in  four  sections  appears  to  range  from  1-4 
percent  and  total  pore  porosity  for  all  sections  ranges  from 
less  than  a  percent  to  roughly  20  percent .  Spheroidal  pores 
dominate;  1)  evenly  dispersed  limonite-f i 1 led  vugs  the  size 
of  f oraminif era ,  and  2)  equidimens iona 1  pores  with  indistinct 
boundaries  barely  visible  to  the  unaided  eye.  Other  than 
limonite  in  concentrations  less  than  0.5  percent,  no  sec¬ 
ondary  minerals,  such  as  silica  or  clay,  were  observed  in 
crack  sections.  Therefore,  based  on  the  theory  predictions 
of  Toksdz  et  al.  (1976)  and  crack  section  analysis,  existing 
microfractures  within  the  Santa  Elena  rock  matrix  must  be 
healed . 
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In  conclusion,  laboratory  results  corroborate  1)  the 
observations  made  from  the  hand  specimens  (with  the  unaided 
eye)  that  the  Santa  Elena  is  composed  primarily  of  calcite 
and  has  a  variable  connected  porosity  range  of  0^  to  10 
percent,  and  2)  indicate  that  although  connected  porosity 
varies  from  2  to  14  percent  in  the  hand  specimens,  the  Santa 
Elena  rock  matrix  is  homogeneous  at  the  macroscopic  scale. 
Randomly  distributed  pore  space,  ellipsoidal  to  spheroidal  in 
geometry  with  radii  less  than  0.5  mm,  and  absence  of  open 
air-filled  fractures  at  the  microscopic  scale  in  the  Santa 
Elena,  play  a  role  in  slightly  reducing  matrix  P-  and  S-wave 
velocities.  Clay  is  also  a  contributing  factor  to  the  slight 
velocity  reduction. 

The  in  situ  seismic  data  of  Reinke  and  Logan  (1983)  and 
Golden  et  al.  (1985)  cannot  be  explained  by  these  models 
because  the  laboratory  P-wave  velocity  measurements  exceed 
the  in  situ  velocity  measurements  by  approximately  40  to  70 
percent.  Our  next  task  is  to  determine  the  cause  this  veloc¬ 
ity  difference. 
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V.  Application  of  Fracture  Model  Theory 


The  main  objective  of  this  research  is  to  investigate 
the  effects  macrofractures  have  on  in  situ  seismic  wave 
velocities.  In  this  chapter,  we  examine  the  relationship 
between  several  measurements  of  macrofracture  orientation, 
porosity,  and  density  and  in  situ  P-,  Sv-,  and  Sh-wave  veloc¬ 
ities  at  varying  azimuths.  Matrix  and  whole  rock  parameters, 
estimated  from  the  laboratory  measurements  of  P-  and  S-wave 
velocities,  pore  porosity,  dry  bulk  and  saturated  rock  densi¬ 
ties,  and  grain  density  for  the  Santa  Elena  limestone,  aid  in 
examining  this  relationship.  Laboratory  and  in  situ  velocity 
measurements  provide  the  parameters  necessary  for  the  appli¬ 
cation  of  theoretical  models  relating  pore  porosity,  fracture 
porosity,  and  fracture  density,  to  seismic  wave  velocities. 

A  discussion  of  these  various  theoretical  models  immedi¬ 
ately  follows.  Then,  model  parameters,  determined  fromi  in 
situ  velocity  and  surface  fracture  measurements  and  labora¬ 
tory  data  (from  chapter  IV),  are  given.  Next,  application  of 
Thomsen's  "Biot-Consistent"  model  (1985),  using  solid  grain 
and  Biot  medium  parameters  and  in  situ  P-  and  S-wave  velocity 
measurements,  predicts  total  porosity  and  fracture  density 
along  several  seismic  refraction  lines.  These  predictions 
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are  then  compared  to  the  in  situ  model  parameters  determined 
from  surface  fracture  measurements. 

Model  Theory 

Several  theories,  such  as  the  "Noninteracting" 

(Budiansky  and  O'Connell,  1980;  Kuster  and  Toksoz,  1974), 
"Augmented  Self-Consistent"  (O'Connell  and  Budiansky,  1977) 
and  "Biot-Consistent"  (Thomsen,  1985) ,  model  the  elastic  mod¬ 
uli  of  porous  rocks. 

According  to  Thomsen  (1985) ,  these  theories  share  the 
minimal  assumptions  concerning  the  structure  of  pore  space  in 
Biot-Gassmann  theory.  Biot  (1962)  derived  the  basic 
equations  for  a  porous,  linearly  elastic,  isotropic  aggregate 
at  low  frequency.  Assuming  that  pore  space  is  interconnected 
and  that  the  frequency  is  sufficiently  low  such  that  the  pore 
fluid  pressure  is  uniform  within  a  given  unit  volume  V,  these 
constituent  equations  may  be  written  as 

(5.1a)  X  =  [I  y, 

(5.1b)  p  =  -K  0  +  ap  ,  , 

p  f  =  ~aM  Q  +  ^  . 

(5.1c)  /V 
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The  bars  indicate  a  volumetric  change  over  V  ,  fixed  in  a 
solid  framework  containing  many  grains  and  pores.  The  shear 
stress  and  pressure  p  are  averaged  entirely  over  V  for  both 
solid  and  fluid  portions.  The  shear  strain  and  dilitation  ^ 
are  also  averaged-  over  V  .  The  incremental  fluid  pressure 
and  incremental  pore  volume  are  pf  and  AVp  respectively;  and 
elastic  parameters  [1*,  K* ,  a,  and  M  are  functions  of  the 
stresses  and  strains  and  fluid  pressure  of  the  initial  state. 
Equations  (5.1)  implicitly  assume  that  ;he  solid  parts  of  '/ 
are  homogeneous  and  isotropic  on  a  microscopic  scale. 

Gassmann  (1951a)  interprets,  without  derivation,  the 
elastic  parameters  in  terms  of  the  solid  and  the  pore  space 
separately.  Because  Gassmann  .made  a  nontrivial  extension  of 
Biot's  work  for  elastic  parameters,  the  r<=>suli-S  are  referred 
to  as  the  Biot-Gassmann  formulas .  When  the  medium  is  in  a 
drained  state,  equations  (5.1a)  and  (5.1b)  reduce  to  Hooke's 
equations  with  shear  and  bulk  elastic  moduli  (or  "frame 
moduli")  ,  |l*  and  K*,  respectively. 

For  connected  pore  space,  Geertsmia  (1957)  and  Mur  and 
Byerlee  (1971)  show  that 
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where  K5  and  <])  are  the  incompressibility  of  the  solid  grains 
and  total  porosity,  respectively. 

The  Biot-Gassmann  formulation  (1941)  makes  two  well- 
known  predictions  for  elastic  moduli  of  porous  rocks:  1)  the 
shear  modulus  of  an  unsaturated  rock  (permeated  by  a  com¬ 
pressible  fluid,  e.g.,  gas)  equals  that  of  the  same  rock 
saturated  with  liquid  and  2)  the  unsaturated  and  saturated 
bulk  modulus  differ  by  a  defined  amount.  The  theory  of 
Toksoz  et  al.  (1974),  developed  for  small  values  of  porosity 
and  fracture  density  ([c(am)]/Otm  <  1),  is  not  consistent  with 

the  latter  of  these  two  predictions. 

Similar  to  Toksoz  et  al .  (1974)  but  in  agreement  with 
both  Biot-Gassmann  predictions,  the  noninteracting  theory  of 
Budiansky  and  O'Connell  (1980)  combines  the  standard  theory 
for  the  elasticity  of  a  solid  isotropic  matrix  with  a  dilute 
concentration  of  spherical  pores  (Eshelby,  1957)  and  the  cor¬ 
responding  theory  for  dilute  concentrations  of  thin,  ellip¬ 
soidal  fractures  (Bristow,  1960) .  The  Budiansky  and 
O'Connell  model  (1980)  assumes  1)  fluid  pressure  equalization 
between  the  two  populations,  and  2)  a  surrounding  medium 
identical  to  that  of  the  solid  grains.  Thus  no  elastic 
interaction  takes  place  between  neighboring  heterogeneities 
spaced  far  apart.  When  no  elastic  interaction  takes  place 
the  effects  on  the  moduli  are  additive  and  lead  to  relation¬ 
ships  linear  i.n  porosity  or  in  fracture  density.  However, 
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because  high  values  of  porosity  and  fracture  density  violate 
its  assumptions,  the  theory  is  not  applicable  to  the  Lajitas 
data.  (The  next  section.  Model  Parameters,  gives  fracture 
porosity  estimates  and  fracture  density  estimates  determined 
through  the  application  of  equation  (5.41)  by  using  fracture 
length  and  width  parameters  measured  at  the  surface  in 
outcrop . ) 

The  augmented  self-consistent  model  (cf.,  Budiansky  and 
O'Connell,  1976;  O'Connell  and  Budiansky,  1974,  1977;  Berry¬ 

man,  1980)  statistically  calculates  the  interaction  between 
neighboring  inhomogeneities.  This  theory  assumes  that  the 
solution  of  the  noninteracting  Budiansky  and  O'Connell  (1980) 
model  for  a  single  pore  enveloped,  not  by  the  solid,  but  by  a 
uniform  medium  with  the  elastic  properties  of  the  "whole 
rock"  yields  the  effect  of  many  spherical  pores.  This  model 
allows  for  large  porosity  and/or  fracture  density  but  only 
agrees  with  Biot-Gassmann  theory  in  the  case  of  no  fractures. 

Thomsen  (1985)  provides  the  mathematical  relationships 
and  discusses  the  Biot-Gassmann  theory  and  noninteracting  and 
augmented  self-consistent  models  in  greater  detail.  He  also 
proposes  a  "Biot-Consistent"  model,  applicable  at  low  fre¬ 
quencies  with  no  limits  on  fracture  density  and  total 
porosity . 

The  Biot -Cons istent  model  includes  a  third  dependent 
state  variable  not  included  in  the  non  interact ing  and  self- 
consistent  models  but  recognized  in  Biot-Gassmann  theory--the 
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fluid  pressure  pf.  In  the  undrained  case  pf  is  equivalent  to 
the  incremental  pore  volume  AVp  /V  .  (See  equations  (5.1)  .) 

The  corresponding  modulus  can  be  defined  as  the  pore  incom¬ 
pressibility  Kp  : 


(5.3) 


K' 


Pf 


Consequently,  three  rather  than  two  characteristic  moduli 
(}l,  K,  Kp)  specified  at  any  particular  value  of  K'  (the  satu¬ 
rating  fluid  bulk  modulus)  define  the  elastic  response  at  all 
saturations  for  a  "Biot  medium".  Hence  in  the  model  theory 
for  a  porous  rock,  the  surrounding  medium  of  the  noninteract¬ 
ing  model  possesses  the  three  moduli  of  the  Biot  medium. 

Considering  a  Biot  medium  with  a  substantial  equant 
porosity  and  fracture  density,  all  interconnected,  Thomsen 
(1985)  relates  the  shear  and  bulk  moduli  H  and  K,  respec¬ 
tively,  to: 
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Defining  the  above  parameters  in  terms  of  Poison's  ratio  of 
the  Biot  medium  (subscript  B)  : 

a  -  ^ 

as  -  —I - r  , 

(5.4c)  ~  ''b) 

and 


(5 . 4d) 


for  the  pores,  and 


(5.4e) 

and 
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for  the  fractures.  These  functions  depend  on  Poison's  ratio 
Vb  for  the  Biot  medium: 


(5.4g) 


For  the  shear  and  bul)c  modulus  of  the  Biot  medium,  pb  and  Kb 


(5.4h) 


=  (vs/  Pb  . 


and 


(5.4i) 
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where 


(5.4j)  Pb  ^JPs 

and  (t)B  is  the  porosity  of  the  Biot  medium  and  Pf  and  ps  are 

the  density  of  the  fluid  and  solid  grains  respectively.  The 
total  porosity  4)  equals  the  sum  of  the  total  fracture 
porosity  <t)c  and  pore  porosity  4)p 

(5.4k)  = 


and  the  fracture  density  E  is  related  to  (t)c  (for  circular, 
i.e.,  penny-shaped  fractures)  by 


{5.41) 
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where  the  thickness/diameter  of  the  fractures  equal  the 
aspect  ratio  X.  For  a  spectrum  of  fracture  shapes  the  aver¬ 
age  fracture  density  over  the  spectrum,  <£>,  replaces  E  in 

equations  (5.4a)  and  (5.4b)  .  The  shear  and  bulk  moduli  of 
the  solid  grains,  P-s  and  Ks  ,  respectively  equal: 
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(5 . 4m) 


P3, 


and 

Ks  =  (VP,)'  p,  -|ivs/  p^. 

(5 . 4n)  ^ 

In  the  drained  case  (pf  =0)  ,  as  according  to  the  Biot-Gassmann 
theory  that  the  shear  modulus  drained  equals  the  shear  modu¬ 
lus  saturated 

(5.5a)  ' 

and  the  in  situ  bulk  modulus  K(pf  )  reduces  to 


K(pf  )  =  K(0)  =  K,  1  - 


(5.5b) 


1  -  d» 


-  A*b  <1)  =  K 


where  "starred"  parameters,  such  as  a*B/  refer  to 

the  Biot  medium  in  the  drained  case. 

Note  that  in  the  drained  state,  the  equations  of  motion 
define  elastic-wave  velocities  in  terms  of  moduli  and  density 
p*: 


(5.6a) 


and 
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vs*  = 


(5.6b) 


with  the  densitv  given  by 


(5.6c) 


P*  =  0p,+  (l-(f 


(The  quantities  VP*  and  VS*  represent  in  situ  seismic  wave 
velocities.)  Substituting  equation  (5.6c)  into  equations 
(5.6a)  and  (5.6b)  and  solving  for  )i*  and  K*  : 


(5.7a: 


p*  =  (vs*)V,  -  (vs*)V,<l), 


(5.7b) 


K*  =  (VP*)p^  -  ^VS*)p^  +  ^VS*rp,  -  (VP*)P3  0 


Equating  equations  (5.7b)  to  corresponding  equation  (5.5b) 
and  simplifying: 
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Thus,  application  of  the  Biot -Consistent  model  using  parame¬ 
ters  determined  from  geophysical  field  measurements  (VP*  and 
VS*)  and  laboratory  measurements  (ps.  Kg,  Ps,  a*B,  b*B,  A*b, 

B  B/  arid  ({ip)  predicts  fracture  density  and  fracture  porosity 
(0  =  0p  +  <!>c)  iri  the  field. 

Model  Parameters — In  Situ  Measurements 

The  main  objectives  of  the  geophysical  field  work  were 
1)  to  measure  in  situ  P-  and  S-wave  velocities  at  varying 
azimuths  using  shallow  refraction  seismic  methods  and  2)  to 
select  unit  volumes  at  or  near  the  chosen  refraction  site  and 
count  and  measure  the  length  and  orientation  of  each  observed 
fracture  trace.  Fracture  density  and  fracture  porosity  esti¬ 
mates  made  from  the  surface  fracture  length  and  width 
measurements  give  a  comparison  to  those  predicted  from  the 
Biot-Consistent  model  using  in  situ  velocity  measurements  and 
parameters  derived  from  the  laboratory  data  for  the  Santa 
Elena  limestone.  Determining  P-  and  S-wave  velocities  along 
refraction  spreads  aligned  parallel  and  perpendicular  to  the 
strikes  of  the  dominant  observed  fracture  planes  tests  for 
azimuthal  anisotropy. 

Data  Acquisition.  In  this  study,  an  iron  plate  and  ham¬ 
mer  provided  both  the  compressional  and  shear  source.  The 
compressional  source,  according  to  the  representative  sam¬ 
pling  of  field  seismograms  shown  in  Appendix  D,  generated  a 
range  of  frequencies  from  200  -  300  Hz  and  wavelengths  from  9 


73 


74 


-  15  m  in  the  Santa  Elena  limestone.  The  shear  source 
(Figure  5.1)  used  in  these  experiments  transmitted  a  maximum 
horizontal  range  of  approximately  50  m  (as  opposed  to  a  100  m 
horizontal  range  for  the  compressional  source)  and  generateri 
frequencies  from  75  -  330  Hz  and  wavelengths  from  5  -  35  m. 

To  observe  slight  anisotropic  effects  in  fractured 
isotropic  and  homogeneous  material,  the  source  must  generate 
wavelengths  much  greater  than  the  fracture  size  and  separa¬ 
tion  distance  between  fractures  (Backus,  1962)  .  Wavelengths 
ranging  from  5  to  35  m  are  not  much  greater  than  the  35  - 
40  m  separation  distance  measured  between  macrofractures  on 
the  lineation  map.  However,  the  2.5  -  30.5  cm  range  of 
separation  distances  measured  between  fractures  observed  in 
outcrop  ranges  from  0.0007  -  0.007  times  the  size  of  these 
wavelengths.  These  smaller  scale  fractures  (overlapping  the 
lower  macrofracture  and  upper  microfracture  boundary-  scale 
limits)  vary  0.1  -  0.3  cm  in  width  and  1  -  110  cm  in  length. 
Consequently,  the  high  frequency  content  and  short  range  of 
the  shear  and  compressional  source  constrained  the  study  to 
the  examination  of  the  effects  that  the  fractures  observed  in 
outcrop  have  on  in  situ  P-  and  S-wave  velocity. 

Figure  5.2  shows  the  geometrical  arrangement  of  refrac¬ 
tion  lines  used  to  measure  P-  and  S-wave  velocities  at  sev¬ 
eral  azimuths  and  locates  the  unit  volumes  chosen  for  the 
fracture  orientation  and  density  and  fracture  porosity 
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F igu  re  5.2  Map 
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determinations.  The  refraction  lines  arranged  in  a  "star- 
like"  fashion  aided  in  distinguishing  azimuthal  velocity 
variations  due  to  fractures,  from  velocity  variations  due  to 
stratigraphic  discontinuities  (Telford  et  al  .  ,  1,976;  Crampin 
1984b) ;  whereas  reversed  refraction  spreads  facilitated  the 
detection  of  shingling  effects  already  observed  at  the  site 
(Reinke  et  al . ,  1983)  .  Minimizing  effects  of  shear-wave 
splitting  by  orienting  refraction  lines  with  major  fracttre 
trends  helped  in  recognizing  the  effects  of  anisotropy 
(Crampin,  1978,  1981,  1984a,  1985,  Crampin  et  al.,  1984).  In 
addition,  measurement  of  P-  and  S-wave  velocities  perpendicu¬ 
lar  and  parallel  to  dominant  air-filled  fracture  sets  tested 
for  azimuthal  anisotropy  (Lynn  and  Thomsen,  1986;  Rai  and 
Hanson ,  1986). 

Fracture  orientations  measured  in  unit  volumes  A  -  E 
trend  to  the  NNW  and  ENE  (Figure  5.3a-e)  .  Figure  5.3f  shows 
the  relationship  between  the  six  refraction  line  orientations 
and  the  cumulati^/e  representation  of  fracture  trace  orienta¬ 
tions  observed  in  unit  volumes  A  -  E.  Excluding  spread  4, 
all  refraction  lines  trend  subparallel  to  one  of  the  observed 
sets  of  fracture  planes. 

For  spread  1,  high  resolution  profiles  were  obtained  by 
using  a  2  m  geophone  spacing  and  a  12  channel  ana  log-to-dig- 
ital  recorder  (sampling  at  10  |ls  intervals),  and  by  making 

several  in  line  shots  offset  appro.ximately  22  m  from  one 
anot.her  on  both  sides  cf  the  spread.  The  mobile  source  and 


77 


Figure  5.3a-e  Observed  strikes  of  fracture  planes  in 
unit  volumes  a)  A,  b)  B,  c)  C,  d)  D  and  e)  E. 
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Figure  5.3f  Illustration  showing  refraction  line 
orientation  versus  a  cumulative  representation  of 
observed  fracture  orientations  in  unit  volumes  A  -  E. 
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stationary  receivers  allowed  for  substantial  horizontal  cov¬ 
erage  in  a  short  period  of  time.  The  multiple  shotpoints  for 
each  spread  also  aided  in  resolving  the  ambiguities  inherent 
in  the  seismic  refraction  travel  time  curves  based  on  first 
arrivals  (Ackermann  et  al.,1986) .  Limited  time  and  take-out 
cord  between  channel  hook-ups  constrained  the  radial  array, 
comprised  of  sp''eads  2,  3,  4,  5  and  6,  to  two  geophones  per 
arm.  In  each  radial  arm,  spacing  be  ^en  geophones  equaled 
6  m.  Multiple  shotpoints,  placed  10,  15,  and  20  m  from  the 
nearest  geophcne  in  each  radial  spread,  yielded  si.x  first 
arrival  picxs  per  line. 

Vertical  10  Hz  PE3  Sensorphones  and  horizontal  14  Hz 
Mark  Product  phones  measured  velocity  along  each  spread.  The 
horizontal  phones  measured  shear  wave  velocities  polarized 
parallel  and  polarized  perpendicular  to  each  S-wave 
propagation  path. 

Rocky  surface  conditions  at  the  site  necessitated  the 
burial  of  horizontal  phones  in  small  holes  filled  with  fast 
quenching  criiit  (Krohn ,  1 984 )  .  Using  nuts  and  bolts  rather 
than  spikes,  horizontal  phones  were  oriented  and  leveled  in 
wet  grout  a;:.;  later  re.moved  by  loose.ning  the  nuts  (Figure 
5.4) .  The  Mark  Product  phones  required  24  pitch,  3/8"  diame¬ 
ter  nuts  and  bolts,  with  bolts  measuri.ng  2-3  inches  in 
length.  The  fast  quenching  grout,  termed  "Jug  Plug,  the 
Mighty  Miracle  Mix",  set  within  15  minutes.  Two  and  three- 
quarter  pounds  of  type  I  or  III  Portland  cement,  0.69  lbs  of 
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Horizontal  phone 


Grout 


Figure  5.4  Illustration  of  method  used  to 
reduce  geophone  ground  coupling  at  lew  f re 
quencies  . 
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CalSeal  (gypsum)  ct  ■'nt,  3.44  lbs  of  clean  masonry  sand,  and 
1.72  lbs  (or  0.21  gals)  of  cool  water  made  1/15  ft^  (3"  X  6" 

X  6")  of  "Jug  Plug''.  Developed  for  the  rapid  construction  of 
seismometer  vaults  and  good  long  term  leveling  of  instru¬ 
ments,  this  grout  has  sufficient  viscosity  to  cause  the  mix¬ 
ture  to  seek  its  own  level  when  poured  (Lewis,  1987) .  Hence, 
the  vertical  phones  were  placed  on  small  grout  pods  using 
triangular  stands  (or  TIF  snow  bases) .  Due  to  limited  field 
time,  geophones  were  only  grouted  for  spread  1.  Note  that 
geophone  spikes  penetrated  the  surface  more  easily  in  the 
region  of  the  radial  array. 

P-  and  S-Wave  Velocity.  Using  refraction  techniques, 
Reinke  and  Logan  (1983)  measured  P-wave  velocity  along  an 
1100  m  line  oriented  subparallel  to  and  possibly  overlapping 
spread  1  (and  5)  of  this  study.  Reinke  and  Logan  (1983) 
spaced  24  geophones  at  33.5  m  intervals  and  used  an  8  lb 
Kine-Stick  two  component  explosive  source,  buried  just  below 
the  surface,  and  offset  33.5  and  300  m  from  the  nearest 
geophones  for  forward  and  reverse  profiles  respectively. 

They  then  applied  a.n  automated  seismic  refraction 
interpretation  progra.m,  (SIPT)  to  determine  the  number  of 
subsurface  layers  represented  by  plotted  travel  times  picked 
to  the  nearest  millisecond.  This  program,  developed  by  the 
US  Bureau  of  Mines,  uses  regression  and  ray  tracing 
techniques  to  proouce  a  subsurface  dept.h-velocity  profile 
(Scott,  19'73)  .  A  jagged  interface  between  average  velocities 
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of  2.907  km/s  and  3.779  km/s  resulted  in  an  attempt  by  SIPT 
to  fit  large  time  gaps  present  in  the  arrival  times  from  the 
forward  shot.  Because  of  the  fairly  simple  near  surface 
stratigraphy  and  structure,  Reinke  and  Logan  (1983)  speculate 
that  "shingling",  as  described  by  Spencer  (1965),  causes  the 
observed  time  gaps.  Shingling  occurs  when  "peaks  and  troughs 
move  forward  through  the  envelope  which  defines  the  refracted 
arrival.  In  this  process,  the  amplitude  of  the  first 
extremum  decreases  and  it  is  eventually  lost  in  the  noise. 

At  this  offset  where  extremum  is  lost,  there  is  a 
discontinuity  in  the  time-distance  curve  and  a  new  shingj.e  is 
added  corresponding  to  a  later,  larger  amplitude  extremum" 
(Spencer,  1965)  .  However,  this  phenomenon  does  not  explain 
the  absence  of  time  gaps  in  their  reversed  profile.  Disre¬ 
garding  the  forward  profile  as  invalid  and  using  only  travel 
times  from,  the  reversed  profile  for  interpretation  yielded  no 
significant  P-wave  velocity  interface  in  the  Santa  Elena 
limestone.  The  P-wave  velocity  averages  approximately  3.477 
km/s,  as  indicated  for  the  Santa  Elena  by  the  reversed  shot. 

Travel  time  curves  representing  P-  and  Sv-wave  first 
arrivals  for  spread  1  of  this  study  behave  similarly  (Figures 
5.5  -  5.7) .  The  forward  and  reversed  profiles  for  spread  1 
are  laid  out  in  the  same  direction  as  the  1100  m  forward  ard 
reversed  profiles  of  Reinke  and  Logan  (1983)  .  "Shingling" 
occurs  in  all  profiles.  Understanding  of  this  shingling 
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SV-WAVE  TRAVELTIME  CURVES 
(Spread  1) 


SH-WAVE  TRAVELTIME  CURVES 
(Spread  1) 


FORWARD  PROFILE 


a 


Discance  (m) 


a) 


REVERSE  PROFILE 


Distance  (ml 


b) 

Figure  5.7a-b  Sh-wave  first  arrival 
picks  for  a)  forward  and  B)  reversed 
profiles  of  spread  1 
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phenomenon  requires  further  study  outside  the  scope  of  this 
research . 

Static  effects,  such  as  those  resulting  from  topographic 
gradients,  may  also  cause  time  gaps.  Source  offsets  along  a 
slight  uphill  gradient  extending  from  SI  to  S7  produced  time 
gaps.  Using  expressions  for  dipping  beds  with  discrete 
velocities  (Dobrin,  1976),  one  can  obtain  the  change  in 
elevation  and  the  velocity  of  the  refracting  medium  from  the 
direct  arrival  (Vq)  and  the  forward  and  reversed  slowness  and 
intercept  times  of  the  refractor  (assuming  the  refractor  is 
horizontal)  .  Using  rough  intercept  time  estimates  (±4  ms)  , 
these  expressions  approximate  an  average  change  in  elevation 
of  1.6m.  A  change  in  elevation  of  1.5  m,  determined  in  the 
field  using  hand  leveling  techniques  with  a  Brunton  compass, 
comes  within  10  percent  of  this  1.6m  estimate.  The  average 
refractor  depth  ranges  from  3.2  ±  1.3  m  (or  10.5  ±  4.2  ft)  at 
shotpoint  S3  to  4.8  ±  1.9  m  (or  15.7  ±  6.2  ft)  at  shotpoint 
S7 .  First  arrival  times  pic)ced  to  the  nearest  ms  yield  P-, 
Sv-,  and  Sh-wave  velocities  equal  to  3.02  ±  0.6  Icm/s,  2.53  ± 
0.5  >:m/s,  and  2.24  ±  0.4  )cm/s  respectively.  Broadened  pea)cs, 
due  to  irregularity  and  weathering  of  the  rock  surface,  and 
variations  in  horizontal  and  vertical  velocity  of  overburden 
may  have  introduced  error  in  the  first  arrival  travel  time 
picks  (Domzalski,  1956)  .  Over  relatively  short  horizontal 
distances  (22  -  44  m>  ,  first  anlval  Limes  picked  from  broad 
peaks  to  the  nearest  millisecond  introduced  as  much  as  20 
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percent  error  in  the  velocity  measurements.  In  contrast, 
travel  times  picked  to  the  nearest  millisecond  and  measured 
over  greater  horizontal  distances  produced  less  than  1 
percent  error  in  the  velocity  measurements  deterijiined  for  the 
reversed  1100  m  line  by  Reinke  and  Logan.  Also,  narrower 
pulse  widths  produced  by  the  explosive  source  made  first 
arrival  picks  more  obvious  in  their  study. 

In  the  radial  array,  the  minimum  and  maximum  compres- 
sional  wave  velocity  measurements  (based  on  error  introduced 
by  first  arrival  picks)  equal  2.53  ±  0.5  km/s  and  3.35  ±  0.7 
km/s;  whereas  those  for  shear  wave  velocity  equal  1.59  ±  0.3 
km/s  and  2.70  ±  0.5  km/s  (Figure  5.8) .  Median  values  equal 
2.99  ±  0.6  km/s  and  2.22  ±  0.4  km/s  for  P-  and  S-wave  veloci¬ 
ties  respectively.  In  situ  P-wave  velocities  are  40  to  60 
percent  less  than  the  median  laboratory  Vp  measurement  of 
5.68  km/s;  whereas  for  S-waves,  the  median  Vs  laboratory  mea¬ 
surement  of  3.10  km/s  and  in  situ  velocities  differ  by  20  to 
50  percent.  Because  a  compositionally  homogeneous  medium 
with  air-filled  pores  (spheroidal  to  ellipsoidal  in  geometry) 
has  relatively  little  effect  on  shear  wave  propagation  (in 
comparison  to  compressional  waves) ,  open  -filled  fractures 

most  likely  cause  the  substantial  drop  in  i  i  P-  and  S-wave 
velocities.  The  Vp/Vs  values  of  1.10  and  1.24  in  line  2  and  3 
seem  low.  However,  an  increase  in  Vp  and  decrease  in  Vs  by  20 
perot-.'.t  r=.ises  these  values  to  1.66  and  1.86.  Ignoring 
error,  Vp/Vs  ranges  from  1.10  to  1.97  for  spreads  2  through  6. 


If  error  is  ignored,  higher  radial  velocities  along 
northeasterly  spreads  (in  comparison  to  their  transverse 
counterparts)  suggest  azimuthal  anisotropy.  Compressional 
waves  propagating  parallel  to  dominant  fracture  trends,  N20W 
and  N55E,  exceed  velocities  measured  along  lines  fringing  the 
range  of  dominant  fracture  trends  (N35W,  N85E>  or  lines  par¬ 
allel  to  no  fractures  at  all  (N35E)  by  approximately  5  to  25 
percent.  P-waves  propagating  to  the  northeast  exceed  those 
propagating  to  the  north  northwest  by  7  to  15  percent  in 
velocity.  Other  work  done  at  the  Lajitas  seismic  station 
suggests  that  incoming  teleseismic  waves  from  the  east  travel 
at  faster  velocities  than  those  coming  from  the  north 
(Golden,  personal  communication) .  Therefore,  easterly  frac¬ 
ture  sets  may  have  slightly  more  effect  on  velocities  than 
northerly  sets. 

However,  the  error  introduced  in  line  1  also  applies  for 
lines  2  through  6.  In  addition,  few  data  points  defining  the 
radial  profiles  (2  points  for  every  shot),  and  minor  static 
problems  originating  from  varying  shotpoint  locations  further 
cloud  the  picture.  High  winds,  a  crude  shear  source,  and 
ungrouted  geophones  in  the  radial  array  contribute  further  to 
the  error.  Thus,  azimuthal  anisotropy  can  not  be  determined 
with  certainty.  At  most,  one  can  infer  that  a  dense  network 
of  open  air-filled  fractures  lower  the  expected  wave 
velocities . 
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Thomsen  speculates  that  a  medium  with  open  air-filled 
orthogonal  fracture  sets,  approximately  equal  in  prominence, 
is  statistically  isotropic  (personal  communication) .  The 
next  two  sections.  Surface  Fracture  Porosity  Estimates  and 
Surface  Fracture  Density  Estimates,  closely  examine  the 
geometrical  and  spatial  distribution  of  fractures  exposed  at 
the  surface  within  the  seismic  refraction  array  boundaries. 

Surto^-e  FLacLure  Porosity  Estimates.  Recall  equation 
(4.24)  : 


N 

n*  1 


Now  assume  a  spectrum  and  the  volume  V,  given  by  surface 
area  A  and  depth  d; 


(5.9) 


4>c=  Z  Z 


n=l 


n=l 


For  vertical  fractures  extending  to  the  unit  volume  depth  d, 
estimating  each  opening  of  fracture  traces  exposed  at  the 
surface  in  terms  of  an  ellipse  with  major  and  minor  axis  an/2 
and  bn/2  yields  a  total  fracture  porosity  0c 


.rr  N  w  N 

71  7t 

0.-= -  2_>  an  t)n  d  =  -  > 

4Ad  n  4A  ^ 


■n  *^n  r 


(5.10) 
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where  and  bn  equal  the  measured  length  and  width  of  a  sin¬ 
gle  fracture  trace  exposed  at  the  surface  (Figure  5.9) .  This 
relationship  holds  when  the  total  fracture  porosity  for  a 
given  unit  surface  area  represents  fracture  poros'ity 
throughout  the  region  of  interest . 

Outcrop  exposure  at  various  locations  within  the  seismic 
refraction  array  boundaries  constrained  coverage  of  each  unit 
surface  area.  The  surface  areas  range  from  2270.96  cm^  to 
13548.36  cm2: 


Unit  A 
Unit  B 
Unit  C 
Unit  D 
Unit  E 


13548.36  cm2 
2270.96  cm2 
5903.21  cm2 
3860.00  cm2 
6350.79  cm2 


Observed  fracture  widths  in  unit  volumes  A  -  E  average 
around  0.2  and  0.3  cm.  Holding  the  fracture  width  constant 
at  0.2  cm,  4)c  equals 


(5.11) 


(0.2}k 

4A 


I- 


0.1571 

A 


r4 

I 


Calculating  the  fracture  porosity  and  aspect  ratio  for 
each  individual  frar-.ture  in  a  given  unit  area  (before  sum¬ 
ming)  yields  a  spectrum  of  fracture  populations; 
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fracture 


Unit  Volume 


Figure  5.9  Illustration  showing  fracture  para¬ 
meters  a,  b,  and  d,  for  a  given  unit  volume. 
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^  '■  n=l,  2,  3,  N 

(5. .12a)  A 

where  for  every  fracture  "n",  ())c(n)  is  associated  with  an 
aspect  ratio  On: 

a,  -  ^  ^  . 

(5.12b)  " 

That  is,  every  fracture  with  length  an  and  width  bn  and 
depth  d,  is  associated  with  a  fracture  volume  anbnd  and  an 
aspect  ratio  bn/an-  For  each  measured  fracture  in  units  A 
through  E,  the  aspect  ratio  and  volume  contribution  were 
obtained  using  equations  (5.12).  The  fractures  range  from, 
0.002  to  0.157  in  aspect  ratio.  For  each  unit  area,  these 
fractures  were  grouped  by  powers  of  ten  in  aspect  ratio 
10~2,  10“^,  and  10*^)  and  their  fracture  volumes  summed  to 
yield  c(am),  the  total  fracture  volume  for  each  group. 

Figure  5.10  associates  the  fracture  porosity  (or  concen¬ 
tration)  c(am)/Ad  with  the  average  aspect  ratio  (bn/an) /N  for 

each  group.  Two  populations  dominate--f racture  populations 
with  average  aspect  ratios  of  0.024  and  0.005  with 
concentrations  ranging  from  0.33  -  1.55  percent  and  1.12  - 
2.87  percent  respectively  (Table  5.1a)  .  The  total  fracture 
porosity  ranges  from  1.88  -  3.93  percent  with  a  median  value 
equal  to  3.05  percent.  Similarly,  at  fracture  widths  equal 
to  0.3  cm,  <{)c  equals 
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Figure  5.10  Concentrat ion  of  crack  populations  having 
crack  widths  equal  to  0.2  cm. 


TABLE  5 . la 

ASPECT  RATIO  VS.  CONCENTRA'^ION  (fracture  width  =  0.2  cm) 


Concentration  (%) 

Aspect  Ratio  Unit  A  Unit  B  Unit  C  Unit  D  Unit  E 


0 . 157 

0 . 00 

0.00 

O.OC 

0.00 

0 . 00 

0.064 

0.03 

0.00 

0 . 00 

0.00 

0.00 

0 . 024 

0 . 64 

1.06 

0.33 

i  .  55 

0.76 

0.005 

2.38 

2^.7 

I.J.Z 

1.79 

L.I2 

3.05 

3.93 

2 . 05 

3 . 34 

1 . 88 

TOTAL 

FRACTURE 

CONCENTRATION : 

MEDIAN 

=  3.05% 
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>'c  ^ 


n-1 


(5.131 


0.2356 


I*. 


n*l 


And  equations  (5.12)  for  bn  =  0.3  cm  become: 


(5 . 14a) 


.(I  0.2356  -  M 

<}),  ln|  =  - -  a^  ;  for  n=l,  2,  3,  ...,  N 


and 


(5 . 14bi 


a  = 

n 


bn  0.3 


‘  n  n 


Calculating  the  aspect  ratio  and  the  volume  of  each  measured 
fracture  using  equations  (5.14)  and  grouping  into  populations 
by  aspect  ratio,  also  yields  two  dominant  fracture  popula¬ 
tions  with  average  aspect  ratios  of  0.024  and  0.005  (Figure 
5.11) .  Concentrations  range  from  0.73  -  3.05  percent  for 
aspect  ratios  averaging  0.024  and  from  0.0  -  3.07  percent  for 
ratios  averaging  0.005  (Table  5.1b) .  With  two  additional 
fracture  populations  present  in  small  qua.ntities  when  frac¬ 
ture  widths  equal  0.3  cm,  total  fracture  porosity  ranges  from 
1.82  -  4.87  percent  with  a  median  value  of  4.59  percent. 
Concentrations  for  fracture  populations  having  aspect  ratios 
of  0.064  and  0.137  range  from  0  -  0.14  percent. 
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Figure  5.11  Concentration  of  crack  populations  having 
crack  widths  equal  to  0.3  cm. 


TABLE  5.1b 

ASPECT  RATIO  VS.  CONCENTRATION  (fracture  width  =  0.3  cm) 


Aspect  Ratio 

Unit  A 

Unit 

Concent  rat  ion 
B  Unit  0 

{%) 

Unit  D 

Unit  E 

0 . 137 

0.02 

0.00 

0.00 

0.00 

0.01 

0.064 

0 . 11 

0.06 

0.04 

0.12 

0 . 14 

0.025 

1.39 

2.31 

0.73 

3.05 

1  .  67 

0.005 

3.07 

2 . 50 

1.^7 

0.00 

4.59 

4.87 

3.08 

4.84 

1.82 

TOTAL 

FRACTURE 

MEDIAN 

CONCENTRATION 
=  4.59% 
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The  relatively  large  concentrations  at  the  aspect  ratio 
of  0.005  violate  the  imposed  restriction  on  the  theoretical 
model  of  Toksoz  et  al.  (1976)  that  c (0^) 70^  <  1.0.  In  all 

cases  but  one,  c  (0 . 005) /O . 005  exceeds  1.0  (Table  5 . 2a-b)  .  In 
addition,  the  Toksoz  model  only  applies  for  low  values  of 
pore  porosity  (less  than  10  percent)  and  dilute 
concentrations  of  fractures  (fracture  densities  less  than 
0.1) .  Note  that  air-filled  partings  along  horizontal  bedding 
planes  also  contribute  to  fracture  porosity.  Because  hand 
specimens  extracted  from  the  Santa  Elena  limestone  average  15 
-  20  cm  in  thickness,  four  to  seven  partings  may  occur  in 
every  meter  with  depth.  Five  horizontal  fractures,  0.3  cm  in 
thickness  and  slicing  through  a  1  m^  volume,  adds  1.2  percent 
to  the  cotal  fracture  porosities  listed  in  Tables  5.1a-b.  In 
Pepper's  Mine,  filled  horizontal  partings  along  bedding 
planes  ranges  from  7  -  30  cm  in  thickness  whereas  fillings  in 
vertical  fractures  only  range  from  0-7  cm.  These  observa¬ 
tions  suggest  that  horizontal  partings  may  be  larger  and  con¬ 
tribute  more  to  fracture  porosity  than  do  the  vertical  frac¬ 
tures  observed  at  the  surface.  For  the  case  where  air-filled 
fracture  porosity  due  to  horizontal  partings  equals  that  due 
to  vertical  fractures,  total  fracture  porosity  ranges  from 
3.64  -  9.74  percent.  Thus  theoretical  models  used  to  predict 
total  porosity  and  fracture  density  in  terms  of  in  situ 
velocities  and  matrix  and  whole  rock  parameters 
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TABLE  5.2a 

ASPECT  RATIO  VS.  C  (On)  /ttn  <  1 
(fracture  width  =  0.2  cm) 


C  (ttn)  /ttn  <  1 

Unit  A 

Unit  B 

Unit  C 

Unit  D 

c(0.157) /0.157 

0.0000 

0.0000 

0.0000 

0 . 0000 

C(0.064) /0.064 

o.oc':7 

0.0000 

0.0000 

0 . 0000 

c(0.024) /0.024 

0.2780 

0.4610 

0.1430 

0.6740 

c(0.005) /0.005 

3.9670 

4.7830 

2.8670 

2.9830 

Unit  E 


0.0000 
0.0000 
0.3300 
1 . 8670 
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TABLE  5.2b 

ASPECT  RATIO  VS.  C  (On)  /Otn  <  1 
(fracture  width  =  0.3  cm) 


C  (ttn)  /ttn  <  1 

Unit  A 

Unit  B 

Unit  C 

Unit  D 

Unit  E 

C(0.157) /0.157 

0.0015 

0.0000 

0 . 0000 

0 . 0000 

0 . 0007 

c(0.064) /0.064 

0.0172 

C.0094 

0.0063 

0.0188 

0.0219 

c(0.025) /0.025 

0.5560 

0.9240 

0.2920 

1.2200 

0.6680 

C (0. 005) /O. 005 

6.1400 

5.0000 

4.6200 

3.3400 

0.0000 
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must  hold  for  large  values  of  fracture  porosity  and  pore 
porosity . 

Surface  Fracture  Density  Estimates.  Fracture  porosity 
(4)c)  was  estimated  for  unit  surface  areas  A  through  E  in  the 

V 

previous  section.  This  section  estimates  the  fracture 
density,  e,  for  each  unit  area. 

For  media  containing  ellipsoidal  fractures  of  various 
aspect  ratio,  Budiansky  and  O'Connell  (1976)  relate  the 
fracture  density  £  to  the  number  of  fractures  per  unit  volume 
Nv  and  the  major  axis,  a,  and  minor  axis,  b,  of  each 
fracture : 


(5.14a) 


where 


(5.14b) 


e  = 


2Nv 

n 


t 


and 


(5.14c) 


A  =  Ttab 


equal  a  fracture  perimeter  and  area  respectively.  In  the 
case  of  circular  fractures  (i.e.,  penny  shaped  fractures), 
equation  (5.14a)  simplifies  to 
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(5.14d) 


e  = 


Lynn  and  Thomsen  (1986)  go  further  and  set  equation  {5.14d) 
equal  to  equation  (5.4i),  the  fracture  density/fracture 
porosity  relationship  used  in  the  Biot -Consistent  theoretical 
model  for  circular  fractures 


where  <a^>  is  the  mean  cube  of  the  fracture  lengths  (or 
traces)  exposed  at  the  surface.  This  expression  relates  more 
closely  to  fracture  permeability  rather  than  fracture 
porosity  (Long,  1983) . 

Both  equations  (5.14d)  and  (5.15)  assume  implicitly  that 
all  fractures  within  a  given  unit  volume  must  be  the  same 
size.  For  this  to  be  true,  as  fracture  length  increases, 
fracture  width  increases  proportionately  such  that  the  aspect 
ratio  X  (=  a=  b/a)  remains  constant.  Because  more  than  one 

fracture  population  exists  in  the  Santa  Elena  limestone,  the 
relationship,  e  =  Ny  <a^>,  is  not  appropriate  for  estimating 

surface  fracture  density  at  the  site.  The  existence  of  more 
than  one  fracture  population  requires  that  the  aspect  ratio 
be  talcen  into  consideration  when  estimating  the  surface 
fracture  density. 
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Using  expression  (5.4i),  the  total  fracture  porosity  4)c 
and  the  mean  aspect  ratio  within  a  given  unit  volume,  X,  can 
be  related  to  an  average  value  of  fracture  density  <£>: 


Using  the  aspect  ratios  and  the  total  fracture 
porosities  determined  for  units  A  through  E  in  the  previous 
section,  equation  (5.16b)  yields  a  mean  aspect  ratio  of  0.063 
for  vertical  fractures  with  6.2  cm  openings.  The  average 
fracture  density  ranges  from  0.07  to  0.15  with  a  median  value 
of  0.11  (Table  5.3).  Similarly,  for  vertical  fractures  with 
0.3  cm  openings,  the  mean  aspect  ratio  equals  0.058;  and  the 
average  fracture  density  ranges  from  0.07  to  0.20  with  a 
median  value  of  0.19.  Consideration  of  horizontal  parting 
raises  the  median  density  values  to  0.23  and  0.38,  assuming 
fracture  porosity  due  to  horizontal  parting  equals  that  due 
to  vertical  fracturing  in  the  region  of  study.  Because 
horizontal  parting  appears  extensive  along  cliff  walls,  its 
contribution  to  fracture  porosity  may  even  be  higher.  The 
fracture  density  estimates  for  the  vertical  fractures  act  as 
a  lower  bound. 
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TABLE  5.3 


AVERAGE  FRACTURE  DENSITY  ESTIMATES  FOR  UNITS  A  -  E 

Average  Fracture  Density  <£> 


Unit 

Fracture  Width  «  0.2  cm 
Mean  Aspect  Ratio  •  0.063 

Fracture  Width  =  0.3  cm 
Mean  Aspect  Ratio  “  0.058 

A 

0.11 

0.19 

B 

0.15 

0.20 

C 

0.08 

0.13 

D 

0.13 

0.20 

E 

0.07 

0.07 
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Most  sedimentary  rocks  have  fracture  densities  greater 
than  0.1;  and  rocks  of  interest  to  the  petroleum  industry 
commonly  have  fracture  densities  of  0.3  or  more  (Thomsen, 
1985) .  Hence,  the  estimated  values  for  fracture  density  at 
the  surface  seem  reasonable. 

Inspection  of  fracture  distribution  suggests  that  the 
ENE  trending  fractures  have  more  prominence  than  the  NNW 
trending.  For  unit  areas  A  through  E,  the  mean  fracture 
lengths  range  from  12.01  -  33.58  cm;  and  the  minimum  and  max¬ 
imum  exposed  fracture  traces  equal  1.27  ±  0.05  cm  and  111.76 
±  0.05  cm  (Table  5.4) .  The  mean  fracture  lengths  for  ENE 
trending  fracture  sets  exceed  NNW  trending  sets  by  approxi¬ 
mately  10  -  25  percent  in  units  B,  C,  and  E,  50  percent  in 
unit  A,  and  equal  one  another  in  unit  D.  In  number,  the  NNW 
and  ENE  fracture  sets  equal  130  and  124  respectively.  The 
average  separation  distance  between  NNW  fractures  is  approxi¬ 
mately  6.65  cm;  for  ENE  fractures  the  separation  distance 
averages  6.15  cm.  Although  the  NNW  and  ENE  fracture  sets 
approximately  equal  one  another  in  number  and  in  average  sep¬ 
aration  distance,  the  relatively  short  NNW  fractures  give  the 
ENE  set  more  prominence.  In  addition,  the  infrared  photo 
reveals  more  pronounced  east  trending  macrofractures  (Figure 
2.3) .  These  observations  give  the  in  situ  velocity  measure¬ 
ments  more  credibility. 
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-shaped  geometry. 


Model  Parameters--Solid  Grain  and  Biot  Medium 


The  moderate  to  high  fracture  porosity  and  fracture  den¬ 
sity  estimates  for  both  the  horizontal  and  vertical  frac¬ 
tures,  and  the  2.5  -  14  percent  range  in  pore  porosity  mea¬ 
sured  in  the  laboratory,  suggest  that  application  of  the 
Biot-Consistent  model  is  most  appropriate  in  analyzing  the 
relationship  between  fracture  density  and  porosity  and  in 
situ  P-  and  S-wave  velocities.  Utilizing  the  Biot-Consistent 
model  to  obtain  predictions  of  fracture  porosity  and  fracture 
density  along  seismic  ray  paths  requires  solid  grain  and  Biot 
medium  parameters  as  well  as  pore  porosity  and  in  situ  P-  and 
S-wave  velocity  measurements. 

When  matrix  velocities  equal  the  P-  and  S-wave  intercept 
velocities  of  6.04  km/s  and  3.23  km/s  in  Figures  4.5  and  4.6, 
and  grain  density  equals  that  of  pure  calcite;  equations 
(5.4m,)  and  (5.4n)  yield  a  solid  grain  shear  and  solid  grain 
bulk  modulus  of  2.838  x  10^^  kg/ (km)  s^  and  6.139  .<  10-^ 
kg/(km)s2  respectively  (Table  5.5). 

Table  5.6  lists  the  Biot  medium  parameters  calculated 
for  various  values  of  pore  porosity  using  equations  (5.4c) 
through  (5.4])  .  All  Biot  medium  para.meters  (the  composite  P- 
and  S-wave  velocities,  shear  and  bulk  moduli,  pore  and  frac- 
tire  param.eters,  etc.)  decrease  as  pore  porosity  increases. 

In  the  next  section,  application  of  the  Bict-Consistent 
model  predicts  fracture  porosity  and  density  using  1)  the 
solid  grai.n  parameters  fis  and  Kg,  2)  the  Biot  .medium  parame¬ 
ters,  A3*,  B3*,  ap*,  and  bg*,  calculated  from  the  mean  and 
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TABLE  5.5 


SOLID  GRAIN  PARAMETER  VALUES 

Solid  Grains  : 

VPg  =  compressional  velocity  of  the  solid  grains 

=  compressional  intercept  velocity  (Figure  4.5) 
=  6.04  km/ s 

VSg  =  shear  velocity  of  the  solid  grains 

=  shear  intercept  velocity  (Figure  4.6) 

=  3.23  km/ s 

Ps  =  density  for  pure  calcite 
=  2.72  X  10^2  kg/km3 

Ps  =  2.833  X  10^3  kg/ (km) s2 

Kg  =  6.139  X  10^3  kg/ (km) s2 
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BIOT  MEDIUM  PARAMETER  VALUES  (pf=0) 
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median  laboratory  measurements  of  pore  porosity  and  P-  and  S- 
wave  velocities,  and  3)  the  in  situ  P-  and  S-wave  velocity 
measurements  VP*  and  VS*. 

Model  Application  and  Interpretation 
Applying  the  Biot-Consistent  model  for  the  drained  case 
with  the  in  situ  velocity  measurements  and  solid  grain  and 
•  Biot  medium  parameters  (discussed  in  the  previous  section) 

yields  a  range  in  average  fracture  porosity  of  10  to  15 
percent  when  pore  porosity  equals  6  percent,  the  mean  of  the 
laboratory  porosity  measurements  (Table  5.7) .  For  the  median 
5.2  percent  pore  porosity  measurement,  the  fracture  porosity 
range  is  slightly  higher.  Fracture  density  ranges  from  0.30 
to  0.48.  These  values  come  well  within  the  range  of  fracture 
density  estimates  made  from  the  surface  fracture  measurements 
of  length  and  width  when  horizontal  fracturing,  in  addition 
to  vertical  fracturing,  is  talcen  into  consideration.  Using 
the  predicted  values  of  average  fracture  porosity  and  density 
in  equation  {5.16a)  predicts  a  range  in  average  aspect  ratio 
,  of  0.06  -  0.07  along  spreads  4  through  6,  and  0.10  -  0.11 

along  spreads  1  through  3.  In  comparison,  the  dverage  aspect 
ratios  used  to  estimate  fracture  density  from  surface 
fracture  measurements  equal  0.058  and  0.060  (Table  5.3). 

In  the  radial  array  (spreads  2  through  6)  the  Biot- 
Consistent  model  predicts  higher  fracture  densities  along  the 
north  trending  lines,  indicating  that  the  ENE  fracture  set  is 
more  prominent.  Along  north  trending  spread  1,  measurement 
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CALCULATED  FRACTURE  POROSITY,  ASPECT  RATIO,  AND  DENSITY 
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of  P-  and  S-wave  velocities  at  a  greater  depth  explains,  in 
comparison  to  the  north  trending  radial  spreads,  the  lower 
fracture  density  prediction. 

Thus  the  Biot-Consistent  model  for  the  drained  case  and 
the  data  collected  in  this  study  can  explain  the  difference 
between  the  bulk  rock  and  in  situ  velocities  at  the  Lajitas 
site  . 
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VI.  Summary  and  Conclusions 


A  breached  fractured  antiform,  caused  by  a  lacolithic 
intrusion  sometime  during  the  Tertiary,  exposes  the  Santa 
Elena  limestone  of  Cretaceous  age.  orthogonal  air-filled 

fracture  sets  existing  in  the  Santa  Elena  strike  NNW  and  ENE 
with  subvertical  to  vertical  inclinations.  Primarily  calcite 
with  secondary  clay  and  silica  in  concentrations  less  than  3 
percent  and  6  percent  respectively,  the  Santa  Elena  rock 
matrix  is  homogeneous  composit ionally .  In  comparison  to  P- 
and  S-wave  velocities  for  pure  calcite  (Vp  =  6.53  km/s,  Vg  = 
3.36  km/s),  a  variation  in  pore  porosity  ranging  from  2.5 
percent  to  14  percent,  with  void  space  spherical  in  geometry, 
lowers  compressional  wave  velocities  10  to  25  percent  and  for 
shear  wave  velocities  5  to  15  percent.  The  P-  and  S-wave 
velocities  for  the  solid  grains,  determined  from  intercept 
velocities,  equal  6.04  km/s  and  3.23  km/s  respectively.  In 
situ  velocity  ranges  from  2.53  km/s  to  3.35  km/s  for  com¬ 
pressional  waves,  and  from  1.59  km/s  to  2.28  km/s  for  shear 
waves . 

The  Biot -Consistent  model  for  the  drained  state  accu¬ 
rately  predicts  fracture  porosity  and  average  fracture  den¬ 
sity  with  the  in  situ  velocity  measurements  and  solid  grain 
and  Biot  medium  parameters  determined  from  the  laboratory 
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measurements  for  the  Santa  Elena  limestone.  Average  fracture 
densities  range  from  0.30  to  0.48;  whereas  crack  porosity 
varies  from  10  to  17  percent. 

Several  observations  suggest  the  presence  of  slight 
azimuthal  anisotropy:  1)  the  faster  in  situ  Sv-wave  velocity 
measurements  of  2.70  ±  0.5  km/s  and  2.28  km/s  ±0.5  km/s 
along  easterly  paths  (as  opposed  to  the  Sh-wave  velocity  mea¬ 
surements  of  2.21  ±  0.4  km/s  and  2.22  ±  0.4  km/s),  2)  greater 
fracture  density  predictions  along  north  trending  seismic  ray 
paths,  3)  longer  more  pronounced  ENE  fracture  sets,  as 
observed  from  the  infrared  photo  and  measured  in  unit  surface 
areas  A  through  E,  and  4)  the  personal  communication  from 
Golden  that  incoming  teleseismic  P-  and  Sv-waves  from  the 
east  travel  at  faster  velocities  than  those  coming  from  the 
north . 

Determining  the  presence  of  azimuthal  anisotropy  at  the 
Lajitas  site  with  certainty  requires  a  refraction  survey 
designed  to  measure  in  situ  P-  and  S-  wave  velocities  with 
less  than  2  percent  error.  Using  an  explosive  source  to  pro¬ 
duce  narrower  pulse  widths,  and  placing  receivers  at  two 
meter  intervals  along  full  spread  lengths  (oriented  at  sev¬ 
eral  azimuths)  to  provide  maximum  horizontal  coverage,  pro¬ 
duces  in  situ  P-  and  S-  wave  velocity  measurements  with 
enough  accuracy  to  determine  the  presence  of  azimuthal 
anisotropy  with  certainty. 
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In  conclusion,  1)  three  sets  of  open  air-filled  frac¬ 
tures  exist  in  the  Santa  Elena  limestone:  horizontal  partings 
along  bedding  planes,  and  vertical  NNW  and  ENE  trending 
fracture  sets,  and  2)  these  open  air-filled  fractures  lower 

V 

the  expected  wave  velocities  at  the  Lajitas  Texas  site. 
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APPENDIX  A 

SCALE  RELATIONSHIPS  AND  ORIENTATION  BETWEEN  MAPS 
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Scale  relationship  between  Figures  1.2,  2.2, 
2.3,  and  3.1. 
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and  location  from  which  the  eight  oriented  han 
specimens  (VLS-1,  VLS-2,  VLS-8)  were  ex¬ 

tracted  . 


APPENDIX  B 

EXPERIMENTAL  SEISMOGRAMS  FOR  SANTA  ELENA  LIMESTONE 
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wave  displacements  polarized  perpendicular  to  horizontal  bedding  planes. 


laboratory  shear  waves  —  PLUG  X1 
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Shear  wave  displacements  polarized  horizontal  and  in  plane  of  ENE  fracture  set. 


UBORATORY  COMPRESSIONAL  WAVES  —  PLUG  Y1B 
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wave  displacements  polarized  horizontal  and  in  plane  of  ENE  fracture  set. 


LABORATORY  SHEAR  WAVES  —  PLUG  Y1B 
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Shear  wave  displacements  polarized  vertical  and  in  plane  of  ENE  fracture  set. 
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wave  displacements  polarized  horizontal  and  in  plane  of  NNW  fracture  set. 
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Shear  wave  displacements  polarized  vertical  and  in  plane  of  NNW  fracture  set. 
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Shear  wave  displacements  polarized  horizontal  and  perpendicular  to  NNW 
fracture  set. 


APPENDIX  C 

SCALARS  Tiijj  AND  'Tijij 
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Defined  below  are  the  scalars  Tajj  and  Tijij  used  by  Toksoz  et 
al..  (1976)  : 


Tijij  -  - 

Fi  ^ 

F3  = 

Fs  ^ 

Fe  =  1  • 


3F- 


2  .  1  ,F4F5  +  F6F7-F8F9 


Tiijj  =  —  +  —  + 


F3  F, 


F2  F, 


where 


=  1  +  A 
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3  ^5.4 

2^  2^~1 
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+  b(3  -  4r)  +  ||a  +  3b)(3  -  4R) 
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+  B0(3  -  4R)  , 


f  A 


1  +  g  -  R(g  +  0)]  +  b(i  -  0)[3  -  4R|  , 
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APPENDIX  D 

FIELD  SEISMOGRAMS  FOR  SANTA  ELENA  LIMESTONE 
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-WAVE  SEISMIC  RECORD  FOR  LINE  1  --  REVERSE 
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University  of  Washington 
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Or.  Stewart  Smith 
IRIS  Inc. 

1616  N.  Fort  Myer  Drive 
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Rondout  Associates 
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Dr.  Jerry  Carter,  Or.  Paul  Pomeroy 
P.O.  Box  224 
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Or.  L.  Sykes 

Lament  Doherty  Geo  logical  Observ. 
Columbia  University 
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Dr.  Pradeep  Taiwan! 

Department  of  Geological  Sciences 
University  of  South  Carolina 
Columbia,  SC  29208 

Or.  R.  B.  Tittmann 

Rockwell  International  Science  Center 

1049  Cam i no  Dos  Rios 

P.O.  Box  1085 

Thousand  Oaks,  CA  91360 

Professor  John  H.  Woodhouse 
Hoffman  Laboratory 
Harvard  University 
20  Oxford  St. 

Cambridge,  MA  02138 

Dr.  Gregory  D.  Young 
ENSCO,  Inc. 

5400  Port  Royal  Road 
Springfield,  VA  22151-2388 


Dr.  Gregory  van  der  Vink 
IRIS 

Suite  1440 
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Dr.  Peter  Basham 
Earth  Physics  Branch 
Geological  Survey  of  Canada 
I  Observatory  Crescent 
Ottawa,  Ontario 
CANADA  K I A  0Y3 

Professor  Ari  Ben-Menahern 
Dept  of  Applied  Mathematics 
Weizman  Institute  of  Science 
Rehovot 
ISRAEL  951729 

Dr.  Eduard  Berg 
Institute  of  Geophysics 
University  of  Hawaii 
Hono lulu,  HI  96822 

Dr.  Michel  Bouchon  -  Uni  vers  ite 

Scientifique  et  Medicale  de  Grenob 
Lab  de  Geophysique  -  Interne  et 
Tectonophys i que  -  I .R. I .G.M-B.P. 

38402  St.  Martin  O'Heres 
Cedex  FRANCE 

Dr.  Hilmar  Bungum/NTNF/NORSAR 
P.O.  Box  51 

Norwegian  Council  of  Science, 

Industry  and  Researcn,  N0R5AR 
N-2007  Kjel ler,  NORWAY 

Dr.  Michel  Campi I lo 
I  .R. I .G.M.-B.P.  68 
38402  St.  Martin  D'Heres 
Cedex,  FRANCE 

C  '.  K  i  n-Y  i  p  Chun 
Geophysics  Division 
Physics  Department 
University  of  Toronto 
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Dr.  Alan  Douglas 
Ministry  of  Defense 
Blacknest,  Brimpton, 

Reading  RG7-4RS 
UNITED  KINGDOM 

Dr.  Manfred  Monger 

Fed.  Inst.  For  Geosciences  &  Nat'l  Res. 

Postfach  510153 

D-3000  Hannover  51 

FEDERAL  REPUBLIC  OF  GERMANY 


Ms.  Eva  Johann isson 
Senior  Research  Officer 
National  Defense  Research  Inst. 

P.O.  Box  27322 
S-102  54  Stockholm 
SWEDEN 

Tor mod  Kvaerna 
NTNF/NORSAR 
P.O.  Box  51 

N-2007  Kjel ler,  NORWAY 

Mr.  Peter  Marshall,  Procurement 
Executive,  Ministry  of  Defense 
Blacknest,  Brimpton, 
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UNITED  KINGDOM  (3  copies) 
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Geophysics  Division 
Geological  Survey  of  Canada 
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NTNF/NORSAR 
P.O.  Box  51 

N-2007  Kjel ler,  NORWAY 
Dr.  Jorg  Sch I i ttenhardt 

Federal  Inst,  for  Geosciences  &  Nat'l  Res. 

Postfach  510153 

D-3000  Hannover  51 

FEDERAL  REPUBLIC  OF  GERMANY 

University  of  Hawaii 
Institute  of  Geophysics 
ATTN:  Dr-  Daniel  Walker 
Honolulu,  HI  96822 
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Or.  Ramon  Cabre,  S.J. 

Observatorio  San  Calixto 
Casi I  la  5939 
La  Paz  Bolivia 

Professor  Peter  Harjes 
Institute  for  Geophysik 
Rhur  University/Bochum 
P.O.  Box  102148,  4630  Bochum  1 
FEDERAL  REPUBLIC  OF  GERMANY 

Or.  E.  Husebye 
NTNF/NORSAR 
P.O.  Box  51 

N-2007  Kjel ler,  NORWAY 

Professor  Brian  L.N.  Kennett 
Research  School  of  Earth  Sciences 
Institute  of  Advanced  Studies 
G.P.O.  Box  4 
Canberra  2601 
AUSTRALIA 

Dr.  B.  Massinon 

Soc i ete  Radi oma na 

27,  Rue  Claude  Bernard 

7,005,  Paris,  FRANCE  (2  copies) 

Dr.  Pi erre  Mech I er 
Soc i ete  Rad iomana 
27,  Rue  Claude  Bernard 
75005,  Paris,  FRANCE 

Dr.  Svein  Mykkeltveit 

NTNF/NORSAR 

P.O.  Box  51 

N-2007  Kjel ler,  NORWAY  (3  copies) 
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Or.  Ralph  A I ew i ne  III 
OARPA/NNRO 

1400  Wilson  Boulevard 
Arlington,  VA  22209-2308 

Or.  Robert  Blandford 
OARPA/NMRO 

1400  Wilson  Boulevard 
Arlington,  VA  22209-2308 

Sandia  National  Laboratory 
ATTN:  Or.  H.  B.  Ourham 
Albuquergue,  NM  87185 

Or.  Jack  Evernden 
USGS-Earthquake  Studies 
345  Middlefield  Road 
Menlo  Park,  CA  94025 

U.S.  Geological  Survey 

ATTN:  Or.  T.  Hanks 

Nat' I  Earthquake  Resch  Center 

345  Middlefield  Road 

Menlo  Park,  CA  94025 

Or.  James  Hanno'i 
Lawrence  Livermore  Nat' I  Lab. 
P.0«  Box  808 
Livermore,  CA  94550 

Paul  Johnson 

ESS-4,  Mail  Stop  J979 

Los  Alamos  National  Laboratory 

Los  Alamos,  NM  87545 

Ms.  Ann  Kerr 
OARPA/NMRO 

1400  Wilson  Boulevard 
Arlington,  VA  22209-2308 

Or.  Max  Koontz 
US  Oept  of  Energy/OP  5 
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&  Engineering 
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U.S.  Geological  Survey 
Ma  i  I  Stop  928 
Reston,  VA  22092 

Or.  Richard  Lewis 
Oir.  Earthquake  Engineering  and 
Geophys i cs 

U.S.  Army  Corps  of  Engineers 
Box  531 

Vicksburg,  MS  39180 

Or.  Robert  Masse' 

Box  25046,  Mail  Stop  967 
Oenver  Federal  Center 
Oenver,  Colorado  80225 

Richard  Morrow 
ACOA/V I 
Room  5741 

320  21st  Street  N.W. 

Washington,  O.C.  20451 

Or.  Keith  K.  Nakanishi 

Lawrence  Livermore  National  Laboratory 

P.O.  Box  808,  L-205 

Livermore,  CA  94550  (2  copies) 

Or.  Carl  Newton 

Los  Alamos  National  Lab. 

P.O.  Box  1663 

Mail  Stop  C335,  Group  ESS-3 
Los  Alamos,  NM  87545 

Or.  Kenneth  H.  Olsen 
Los  Alamos  Scientific  Lab. 

P.O.  Box  1663 

Mail  Stop  C535,  Group  ESS-3 
Los  Alamos,  NM  87545 

Howard  J.  Patton 
Lawrence  Livermore  National 
Laboratory 
P.O.  Box  808,  0-205 
Livermore,  CA  94550 


Mr.  Chris  Paine 

Office  of  Senator  Kennedy 

Sr  3  15 

United  States  Senate 
Washington,  O.C.  20510 

AFOSR/NP 

ATTN:  Colonel  Jerry  J.  Perrizo 
Bldg  410 

Bolling  AFB,  Wash  D.C.  20332-6448 
HQ  AFTAC/TT 

Attn:  Dr.  Frank  F.  Pilotte 
Patrick  AFB,  Florida  32925-6001 

Mr.  Jack  Rach 1 i n 
uses  -  Geology,  Rm  3  C136 
Mail  stop  928  National  Center 
Reston,  VA  22092 

Robert  Reinke 
AFWL/NTESG 

Kirtland  AFB,  NM  87117-6008 

Dr.  Byron  Ristvet 

HQ  DNA,  Nevada  Operations  Office 

Attn:  NVCG 

P.O.  Box  98559 

Las  Vegas,  NV  89193 

HQ  AFTAC/TGR 

Attn:  Dr.  George  H.  Rothe 
Patrick  AFB,  Florida  32925-6001 


Secretary  of  the  Air  Force  (SAFKD) 
Washington,  DC  20330 
Office  of  the  Secretary  Defense 
DOR  &  E 

Washington,  DC  20330 
HQ  ONA 

ATTN:  Technical  Library 
Washington,  DC  20305 

OARPA/RMO/RETR I EVAL 
1400  Wi Ison  Bl vd. 

Arlington,  VA  22209 

DARPA/RMO/Security  Office 
1400  Wi 1  son  B 1 vd. 

Arlington,  VA  22209 

GL/XO 

Hanscom  AFB,  MA  01731-5000 
GL/LW 

Hanscom  AFB,  MA  01731-5000 
DARPA/PM 

1400  Wilson  Boulevard 
Arlington,  VA  22209 

Defense  Techn i ca I 
Information  Center 
Cameron  Station 
Alexandria,  VA  22314 
( 5  cop i es ) 


Donald  L.  Springer 

Lawrence  Livermore  National  Laboratory 
P.O.  Box  808,  L-205 
Livermore,  CA  94550 


Defense  Intelligence  Agency 
Directorate  for  Scientific  & 
Technical  Intelligence 
Washington,  D.C.  20501 


Dr.  Lawrence  Turnbull 
OSWR/NED 

Centra  I  Intelligence  Agency 
CIA,  Room  5G48 
Washington,  D.C.  20505 

Dr.  Thomas  Weaver 

Los  Alamos  National  Laboratory 

P.O.  Box  1663 

MS  C  355 

Los  Alamos,  NM  87545 


Defense  Nuclear  Agency/SPSS 
ATTN:  Or.  Michael  Shore 
6801  Telegraph  Road 

Alexandria,  VA  22310  * 

AFTAC/CA  (ST  INFO)  ; 

Patrick  AFB,  FL  32925-6001 

TACT EC 

Battel le  Memorial  Institute 
505  King  Avenue 

Columbus,  OH  43201  (Final  report  only) 


GL/SULL 

Research  L i brary 

Hanscom  AFB,  MA  01731-5000  (2  copies) 


Mr.  Alfred  I  ieberman 

ACDA/V 1 -OA' State  Department  Building 
Room  5726 

320  -  21St  Street,  NW 
Washington,  D.C.  20451 
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